Special  Programs  for  Analysis  of  Radiation  by  Wire 

Antennas, 


SYRACUSE  UNIV  NY 

01  JUN  1973 


Distribution  Statement  A: 

Approved  for  public  release.  Distribution  is  unlimited. 


Security  Classification 


DOCUMENT  CONTROL  DATA  •  R  &  D 

(Security  daieltlcation  ot  tltlo  ho dy  of  abciract  and  Indtatnf  annofafJut  muif  be  amarad  whan  the  o veratl  upon  It  r/aiif//4dj 

ORIGINATING  ACTIVITY  (Cetporate  cuthot) 

Syracuse  University 

Department  of  Electrical  and  Computer  Engineering 
j  Syracuse,  New  York  13210 

2 •.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

26.  GROUP 

3  REPORT  TITLE 

SPECIAL  PROGRAMS  FOR  ANALYSIS  OF  RADIATION  BY  WIRE  ANTENNAS 

4  oc SCRIPT] VC  NOTH  (Type  ot  report  and  Inclua/va  datee) 

Scientific  Interim 

•  aw  tmoniS)  (Ftret  /*•«*«,  middle  ffljffc.,  laet  (  ‘jme) 

Bradley  J.  Strait 

Tapan  Sarkar 

Dah-Cheng  Kuo 

0  REPORT  OATC 

1  June,  1973 

?«.  TOTAL  NO.  OR  PACES  76.  NO.  OR  REFS 

76  15 

•«.  CONTRACT  OR  GRANT  NO. 

F19628-73-C-0047 

b  PROJECT,  TASK,  ANO  WORK  UNIT  NO. 

5635-06-01 

c.  30D  ELEMENT 

61102F 

4.  OOO  SUBELEMENT 

M.  ORIGINATOR'S  REPORT  MUMBCR(S) 

Scientific  Report  No  1 

•6.  OTHER  REPORT  NO(S)  (Any  other  nunaere  that  may  be  aeelgned 
thie  report) 

AFCRL-TR-73--0399 

10.  DISTRIBUTION  STATEMENT 

A  -  Approved  for  public  release;  distribution  unlimited. 

M  SUPPLEMENT/  RY  NOTES 

TECH,  OTHER 

12.  SPONSORING  MILITARY  AC  TJ  VI  TY 

Air  Force  Cambridge  Research  Laboratories 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts  01730 

ii.  aTTVSact 


Two  user-oriented  computer  programs  are  presented  and  described.  The  first  is 
suitable  for  handling  efficiently  typical  analysis  and  design  problems  involving 
linear  arrays  of  parallel  thin-wire  antennas.  The  second  is  designed  to  enrble 
efficient  analysis  of  radiation  from  vertical  wire  antennas  over  systems  of  radial 
ground  wires.  Examples  are  given  to  illustrate  various  appli:ations  of  both 
programs.  Special  attention  is  devoted  to  use  of  the  first  program  together  with 
a  standard  optimization  procedure  to  design  linear  arrays  of  wire  elements  with 
unequal  spacings  and/or  unequal  wire  lengths. 


t(L' 


DD  ..1473 


UNCLASSIFIED 

Sicunty  Claatificdion 


UNCLASSIFIED 

Security  CUttilice'.ion 

kcv  mono* 

Antennas 
Wire  intennas 
Matrix  Methods 
Optimization 
Radiation 
Array  Analysis 
Array  Design 
Arrays 

Wire  Junctions 
Computer  Programs 
Linear  Arrays 
Radial  Wire  Systems 


AFCRL-TR-73-0399 


SPECIAL  PROGRAMS  FOR  ANALYSIS  OF  RADIATION 
BY  WIRE  ANTENNAS 


by 

Bradley  J.  Strait 
Tapan  Sarkar 
Dah-Cheng  Kuo 


Department  of  Electrical  and  Computer  Engineering 
Syracuse  University 
Syracuse,  Nev  York  13210 


Contract  No.  F19628-73-C-0047 
Project  No.  5635 
Task  No.  563506 
Work  Unit  No.  56350601 


D.  C 

-wnsEfflLQE 


SCIENTIFIC  REPORT  NO.  1 
1  June  1973 


/ 

rVi 


SEP 


10 


Contract  Monitor:  John  F.  Mcllvenna 
Microwave  Physics  Laboratory 


Approved  for  public  release;  distribution  unlimited. 


Prepared  for 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 


SPECIAL  PROGRAMS  FOR  ANALYSIS  OF  RADIATION 
BY  WIRE  ANTENNAS 

Bradley  J.  Strait,  Tapan  Sarkar,  and  Dah-Cheng  Kuo 
Syracuse  University,  Syracuse,  New  York 


ABSTRACT 

Two  user-oriented  computer  programs  are  presented  ?md  described. 

The  first  is  suitable  for  handling  efficiently  typical  analysis  and 
design  problems  involving  linear  arrays  of  parallel  "hin-wire  antennas. 

The  second  is  designed  to  enable  efficient  analysis  of  radiation  from 
vertical  wire  antennas  over  systems  of  radial  ground  wires.  Examples  are 
given  to  illustrate  various  applications  of  both  programs.  Special  atten¬ 
tion  is  devoted  to  use  of  the  first  program  together  with  a  standard  opti¬ 
mization  procedure  to  design  linear  arrays  of  wire  elements  with  unequal 
spaclngs  and/or  unequal  wire  lengths. 
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1.  INTRODUCTION 

This  report  has  three  objectives.  The  first  is  to  present  and  describe 
a  user-oriented  computer  program  suitable  for  handling  typical  radiation 
problems  involving  linear  arrays  of  parallel  thin-wire  antennas.  This  pro¬ 
gram  assumes  that  the  wires  are  all  cf  the  same  radius,  although  they  may 
be  of  different  lengths.  It  is  also  assumed  the  wires  are  centerfed  with 
feed  points  all  in  the  same  plane  and  that  the  wires  are  lossless  with  no 
externally  applied  loading.  The  program  is  based  on  an  analysis  procedure 
suggested  by  Harrington  and  represents  an  application  of  the  method  of 
moments.  [1,2]  Within  Harrington's  general  formalism  the  subsectional 
piecewise  sinusoidal  functions  suggested  by  Richmond  [3]  a^e  used  for  both 
expansion  and  weighting  resulting  in  a  Galerkin  solution  to  the  analysis 
problem.  The  corresponding  computer  program  presented  here  appears  to  have 
an  advantage  with  respect  to  convergence  over  previous  programs,  [4,5]  that 
were  based  on  other  sets  of  expansion  and  weighting  functions.  This  appears 
especially  true  for  radiation  problems  involving  unloaded  thin  wires  near 
resonant  lengths.  The  program  is  described  in  Chapter  2  of  this  report  and 
is  characterized  by  very  simple  data  input  requirements.  The  program  list¬ 
ing  is  included  in  Appendix  B. 

A  second  cbJ  ctive  of  this  report  is  to  demonstrate  use  of  the  program 
described  above  'gether  with  an  appropriate  optimization  procedure  to  treat 
some  typical  optimization  and  design  problems  of  interest.  These  include  de¬ 
sign  of  unequally  spaced  arrays  to  provide  equal  sidelobe  patterns,  selec¬ 
tion  of  wire  lengths  and  interele^ent  spacings  to  optimize  array  directivity, 
and  reduction  of  pattern  sidelobes  through  adjustments  of  wire  lengths. 
Results  for  several  typical  problems  are  included  in  Chapter  3  along  with  a 
discussion  of  the  optimization  procedure  used. 

The  third  objective  of  this  report  is  to  present  a  modified  version  of 
a  user-oriented  computer  program  that  was  presented  in  an  earlier  Scientific 
Report  [4]  together  with  some  results  of  its  application.  This  program  is 
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anc-'her  example  of  use  of  the  Gaierkin  procedure  in  that  subsectional 
piecewise  linear  (triangle)  functions  are  used  for  both  expansion  and 
weighting.  The  initial  program  was  written  specifically  to  analyze 
radiation  and  scattering  by  complicated  configurations  of  arbitrarily 
bent  and  interconnected  thin  wi  s.  As  yet,  however,  data  have  only 
been  presented  for  certain  •  .e-cross  scattering  problems  [5,6]  and  for 
certain  radiation  problems  involving  rectangular  and  circular  wire  loops, 
double-V  antennas,  and  symmetrical-T  antennas.  [5,7]  The  modified  pro¬ 
gram  presented  here  is  useful  for  handling  certain  situations  involving 
special  symmetry  such  as  vertical  wire  antennas  over  radial  wire  (counter¬ 
poise)  systems.  As  mentioned  earlier,  examples  of  its  use  are  included. 
The  program  is  described  in  Chapter  4,  and  the  listing  is  presented  in 
Appendix  C. 


2.  THE  ANALYSIS  PROGRAM 


2.1  Summary  of  the  Method 

The  details  of  Harrington’s  general  method  of  analysis  as  applied 
to  thin-wire  antennas  are  readily  available  elsewhere.  [1,2]  It  is 
assumed  here  the  wires  are  thin  and  that  current  flows  only  in  the  axial 
direction  of  eacv ,  Current  and  charge  densities  are  approximated  by 
filaments  of  current  and  charge  on  che  wire  axes.  The  boundary  condition 
regarding  the  tangential  component  of  the  electric  field  at  the  wire  sur¬ 
faces  is  satisfied  (approximately)  by  requiring  that  the  axial  component 
vanish  at  the  surface  of  each  wire.  In  the  subsectional  approach  used 
here  each  wire  is  thought  of  as  divided  into  a  number  of  short  segments 
or  subsections  connected  together.  The  integrodifferential  equation 
characterising  the  analysis  problem  is  then  reduced  to  a  matrix  equation 
of  the  form 


[V]  =  [Z][I] 


(1) 
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by  expanding  each  wire  current  in  a  sequence  of  subsectional  expansion 
functions  where  each  is  nonzero  only  over  a  small  portion  of  a  wire; 
that  is,  over  a  small  number  of  subsections.  There  are  many  useful 
choices  of  sets  of  expansion  functions,  and  mer.>  nave  been  investigated 
in  detail  [1,2].  For  the  program  pt  isented  here  the  piecewise  sinusoidal 
functions  suggested  by  Richmond  are  used,  and  these  are  depicted  in  Fig.  1.  [3] 


t 


l 


a 

4 


i 


I 


(a)  Sinusoidal  function 


Fig.  1  ~  Subsectional  bases  and  functional  approximation. 

In  order  to  compute  the  elements  of  the  matrices  in  (1)  by  the  method 
of  moments  it  is  necessary  to  define  a  set  of  testing  or  weighting  functions 
in  addition  to  the  current  expansion  functions  used.  If  the  same  functions 
are  used  for  both  expansion  and  test  '.ng  then  the  result  is  known  as  a  Galerkin 
solution  to  the  analysis  problem.  This  is  the  procedure  used  here. 


With  reference  to  Fig.  1  the  nth  current  expansion  function  can  be  written 
as  (the  wires  of  the  array  are  all  assumed  to  be  z-directed) 
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where  uz  is  the  cbvious  unit  vector  and  k  -  2tt/X.  The  field  corresponding 
to  this  single  current  function  is  given  by  [8] 
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where  Ep  and  .1^  are,  of  course,  the  field  components  normal  to  and 
parallel  to  the  direction  of  the  current  element  respectively.  Tf  1^  is 
the  complex  amplitude  of  the  nth  current  expansion  function,  then  the  total 
current  of  the  array  can  be  written  as 

I(z)  **  l  I  I  (z)  (5) 

u  n  n 
n 

and  the  total  field  is  an  obvious  corresponding  generalization  of  (3)  and 
(4)  using  (5). 

Formulas  for  the  elements  of  the  matrices  in  (1)  can  be  found  easily 
from  the  general  theory.  [1,2, 5, 9]  The  matrix  [I]  is  called  the  current 
matrix  and  is  simply  a  column  matrix  of  dimension  equal  to  the  total  number 
of  current  expansion  functions  used.  Each  element  equals  the  complex 
amplitude  of  the  corresponding  current  expansion  function;  that  is,  the 
nth  element  is  1^.  If  the  total  number  of  expansion  functions  used  is 
denoted  by  M  then  [Z]  in  (1)  is  a  square  matrix  of  dimension  M  known  as 


5 


the  generalized  impedance  matrix.  If  E  is  the  vector  field  produced  by  the 

-V  ® 

mth  expansion  function  I^Cz)  then  a  general  formula  for  the  typical  element 
of  [Z]  denoted  by  Z  is 


Z  -  - 
nm 


/  sin  k(z-z  ,) 

|  n*"X  a 

J  sin  k(z  -z  x)  uz 
Zn-1 


E  dz  + 
m 


3ln  . 
sin  k<Vr\)Uz 


E  rfz 
in 


(6) 


Detailed  formulas  for  the  real  and  imaginary  parts  of  Z  are  obtained  by 

nm 

using  (4)  in  (6)  and  are  given  in  Appendix  A  of  this  report  in  terms  of 
sine  and  cosine  integrals.  The  formulas  are  valid  if  all  segment  lengths 
are  equal  along  any  one  wire. 

The  matrix  [V]  in  (1)  is  a  column  matrix  of  dimension  M,  and  its  ele¬ 
ments  are  related  to  the  excitation  of  the  array.  There  is  an  element  of 
[V]  that  corresponds  to  each  expansion  function.  However,  it  is  assumed  here 
that  excitation  voltages  will  be  applied  only  at  the  center  points  of  the  array 
wires.  Hence,  the  only  nonzero  elements  of  [V]  are  those  corresponding  to  cur¬ 
rent  expansion  functions  that  are  centered  at  the  midpoints  or  feedpoints  of 
the  array  wires.  These  nonzero  elements  are  numerically  equal  to  the  complex 
excitation  voltages  that  are  applied.  This  result  is  derived  directly  from 
the  general  formula  for  [V]  by  using  impulsive  excitation  functions  at  the 
points  described. 

Once  the  elements  of  [Z]  and  [V]  have  been  computed  the  current  distri¬ 
bution  is  found  fiom  (1)  by  simply  inverting  [Z],  Input  impedances  and  field 
patterns  of  interest  are  then  determined  using  straightforward  and  well-known 
matrix  manipulations,  completing  the  analysis  problem. 

2.2  Program  Description 

In  the  program  printout  of  Appendix  B  the  main  program  is  listed  first 
followed  by  the  various  subroutines.  The  example  included  in  the  Appendix  is 
for  analysis  of  a  linear  array  of  eight  half-wave  wires  that  are  half-wavelength 


spaced  and  all  ceiterfed  with  real  unit  voltages.  Required  Input  data 
included  in  the  main  program  are  listed  as  follows: 

a)  The  fifth  statement  is 

AK  «  TP  *  (wire  radius)  (TP  *  2n) 

which  inputs  the  wire  radius  in  wavelengths.  (Recall  that  all  array  wires 
are  assumed  to  have  the  same  radius . )  The  wire  radius  used  in  the  example 
included  in  Appendix  B  is  0.007022A. 

fc)  The  total  number  of  wires  making  up  the  array  is  provided  with  the 
sixth  statement  of  the  main  program.  Of  course,  for  the  example  considered 
this  statement  is  simply  N=*8. 

c)  The  eighth  statement  is  the  first  of  (2N-1)  consecutive  data  state¬ 
ments  providing  information  on  the  geometry  of  the  array.  The  first  N  of 
these  a^.e  typified  by  Y(;J)  ■  (half-length  of  the  jth  wire  in  wavelengths)  *  TP 
where  TP  is  3imply  2ir  as  determined  in  statement  four.  The  half-lengths  of 
the  wires  included  in  the  example  of  Appendix  B  ate  all  0.25A.  The  next  (N-l) 
data  statements  indicate  the  positions  of  the  array  vires  relative  to  the 
first.  An  example  of  these  is 

Y(k)  =  (distance  of  the  kth  wire  from  the  first  in  wavelengths)  *  TP 

Note  that  the  wires  in  the  array  of  the  example  are  all  spaced  one  half¬ 
wavelength  apart. 

d)  The  purpose  of  th  DO  LOOP  formed  by  statements  23-25  of  the  main 
program  is  to  provide  information  concerning  excitation  of  the  array.  As 
mentioned  earlier  it  is  assumed  the  wires  are  all  centerfed  so  that  [V]  is 
a  column  matrix  with  each  nonzero  element  equalling  the  complex  excitation 
in  volts  applied  at  the  midpoint  of  the  corresponding  wire.  In  the  illustra¬ 
tive  example  of  the  Appendix  all  wires  are  centerfed  with  real  unit  voltages. 

e)  Part  of  the  computed  output  is  the  field  pattern  of  the  array  taken 
in  the  principal  H-plane.  This  pattern  is  computed  over  180°  of  the  azimuth 
angle  centered  on  the  broadside  direction.  Values  are  calculated  at  intervals 


determined  by  the  26th  statement  of  the  main  program.  The  statement  MQ*3 
shown  in  the  sample  printout  means  that  intervals  of  three  degrees  are 
desired.  If  it  is  only  necessary  to  compute  the  pattern  at  five  degree 
intervals  of  the  azimuth  angle  then,  of  course,  this  statement  should  read 
MQ-5. 

f)  The  only  other  information  required  is  concerned  with  the  number 
of  current  expansion  functions  that  should  be  used  to  represent  the  current 
for  each  array  wire.  This  decision  is  made  by  way  of  statement  AO  included 
in  Subroutine  ASCTFD.  For  analysis  problems  where  an  accurate  representa¬ 
tion  of  the  current  is  needed  to  determine  correctly  input  impedances  and 
other  quantities  of  interest  the  statement 

MM  *  INT(Y(I)  -  0.571)  +  3 

is  recommended,  where  Y(I)  is  defined  in  paragraph  (c)  above.  The  total 
number  of  expansion  functions  for  the  Ith  wire  is  [2(MM)-1]  so  that  the 
relationship  between  this  number  and  the  wire  length  L  is  approximately 

0  <  L  <  |  .  .  .  2(MM)-1  =  5 

•|  <  L  <  0.82A  .  . .  2 (MM)-l  =  7 

0.82A  <  L  <  1.1AA.  .  .  2 (MM)-l  ■  9 
and  so  on. 

This  seems  to  be  a  reasonable  choice  based  on  the  discussion  of  convergence 
included  in  Section  2.3  of  this  report.  If,  on  the  other  hand,  the  primary 
concern  of  the  user  is  with  pattern  characteristics  as  is  the  case  with  the 
optimization  problems  discussed  in  Chapter  3  then  fewer  current  expansion 
functions  are  needed  for  each  wire  and  the  statement 

MM  -  INT(Y(I)  -  0.5/1)  +  1 

can  be  used.  With  this  choice  the  numbers  of  expansion  functions  used  for 

the  intervals  defined  above  are  1,3,5 .  Obviously  the  program  user  can 

make  either  of  these  choices  or  an  independent  one,  depending  on  the  particu¬ 
lar  kinds  of  results  sought. 


Normal  printed  output  of  the  program  includes  all  input  data,  the 
current  distribution  as  expressed  by  the  real  and  imaginary  parts  of  the 
complex  amplitude  of  each  current  expansion  function,  input  impedances  cor¬ 
responding  to  each  feed  point,  and  the  normalized  electric-field  pattern 
in  the  principal  H-plane  together  with  the  normalization  constant. 

Program  operation  proceeds  roughly  as  follows:  The  generalized 
impedance  matrix  [Z]  is  computed  using  instructions  53-97  with  the  obvious 
symmetry  about  the  feed  points  taken  into  account  in  order  to  reduce  overall 
computation  time  and  storage  requirements.  This  matrix  is  inverted  jsing 
Subroutine  CSMIN,  a  complex  inversion  routine  authored  at  the  University  of 
Illinois.  The  current  matrix  [I]  is  determined  by  forming  the  matrix  product 

[I]  -  [Z]"1  [V]  (7) 

in  Subroutine  MULTPY.  Finally,  the  electric-field  pattern  is  computed  for 
the  principal  H-plane  using  instructions  101-120,  and  the  input  impedances 
corresponding  to  feed  points  are  calculated  using  the  DO  LOOP  of  statements 
146-152.  The  FUNCTIONS  ZMN  and  ZM1  compute  mutual  terms  of  [Z]  when  the 
segment  lengths  of  the  corresponding  wires  are  unequal  and  equal  respectively. 
This  is  done  with  the  aid  of  Subroutine  SICI  made  available  through  the  IBM 
Scientific  Subroutine  Package  for  the  purpose  of  computing  the  required  sine 
and  cosine  integrals  S^(x)  and  C^(x).  The  program  described  here  was  developed 
with  an  IBM  370/155  computer  and  storage-time  requirements  are  quoted  with  the 
various  examples  as  they  are  presented  in  later  sections  of  this  report. 

2 . 3  Comments  on  the  Number  of  Expansion  Functions 

The  number  of  current  expansion  functions  that  should  be  used  in  any 
given  situation  has  been  an  open  subject  of  discussion  for  some  time.  It 
is  clear  that  requirements  vary  considerably  in  this  respect  depending  on 
the  kinds  of  results  sought.  For  example,  for  some  pattern  synthesis  and 
optimization  problems  it  may  be  possible  to  use  as  few  as  say  five  expansion 
functions  per  wavelength,  while  In  some  other  situations  where  careful  analyses 
of  current  distributions  are  required  it  is  advisable  to  use  considerably  more, 


say  10-15  expansion  functions  per  wavelength.  Thi'',  of  course,  is  the 
reason  for  incorporating  the  special  instruction  discussed  in  Section 
2f  above. 

Figures  2-5  illustrate  the  kinds  of  results  that  are  obtained  as  the 
numbers  of  expansion  functions  are  varied  for  certain,  sing-e-wire  problems. 
First,  consider  a  half-wave  centerfed  wire  of  radius  0.007022A.  The  input 
admittance  is  calculated  for  this  problem  using  different  numbers  of  current 
expansion  functions  and  for  two  different  choices  of  subsectional  expansion 
functions.  One  of  these  choices  involves  the  piecewise  sinusoidal  functions 
used  here  while  the  other  involves  the  piecewise  linear  (triangle)  functions 
used  for  the  analysis  programs  developed  earlier  to  handle  arbitrary  con¬ 
figurations  of  bent  and  interconnected  wires.  [4]  Computed  results  are  shown 
in  Fig.  2  along  with  the  corresponding  experimental  data  reported  by  Mack.  [10] 
Similar  results  are  shown  in  Figs.  3-5  for  longer  wires.  These  and  other 
available  data  indicate  that  the  piecewise  sinusoidal  functions  offer  somewhat 
faster  convergence  than  the  functions  used  with  the  arbitrary  wire  program  and 
also  that  at  least  9  or  10  expansion  functions  are  needed  per  wavelength  to 
describe  accurately  the  current  distributions  encountered  in  the  kinds  of 
problems  considered  here.  As  indicated  by  Fig.  2  an  even  larger  number  may 
be  advisable  for  wires  near  resonant  lengths. 


Figure  6  provides  computed  information  on  the  self  and  mutual  admittances 
of  two,  half-wavelength,  centerfed  parallel  wires  that  are  spaced  various 
distances  apart  and  are  both  of  radius  0.007022A.  As  before,  data  are  recorded 


as  functions  of  the  number  of  expansion  functions  used  for  each  wire  and  also 
versus  the  separation  distance  d,  for  both  the  piecewise  sinusoidal  and  the 
piecewise  linear  cases.  Once  again,  Mack’s  corresponding  experimental  data  are 
included.  As  before,  the  data  indicate  that  the  piecewise  sinusoidal  functions 
offer  relatively  rapid  convergence  and  that  approximately  10  expansion  functions 
are  needed  per  wavelength  of  wire  to  obtain  good  results  for  the  current  dis¬ 
tributions  and  input  admittances. 
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Fig.  4  -  Input  admittance  of  a  centerfed  linear  antenna  of  length  X 
and  radius  0.007022X  plotted  vs.  n,  the  number  of  sub¬ 
sectional  expansion  functions  used  for  wavelength. 
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The  somewhat  faster  convergence  that  characterizes  the  piecewise 
sinusoidal  current  approximation  for  problems  involving  thin  wires  is 
due  in  part  to  the  inherent  shape  of  the  individual  expansion  functions 
used  and  also  in  part  to  the  fact  that  the  integrations  required  to  ob¬ 
tain  the  elements  of  [Z]  can  be  performed  analytically.  This  is  in  con¬ 
trast  to  certain  programs  based  on  the  piecewise  linear  current  approxi¬ 
mation  that  incorporate  a  substantial  number  of  approximate  numerical 
integrations.  The  piecewise  linear  data  of  Figs.  2-6  were  derived  using 
one  of  these.  It  is  interesting  to  note  that  convergence  is  significantly 
improved  even  in  the  piecewise  linear  case  if  the  integrations  required 
to  obtain  [Z]  are  performed  analytically  [11]. 

Another  advantage  of  the  formulation  used  here  is  that  there  is  no 
tendency  towards  matrix  instability  as  the  number  of  expansion  functions 
is  increased,  a  difficulty  encountered  with  some  other  methods  available 
for  handling  thin-wire  problems. 

3.  OPTIMIZATION  PROBLEMS 

3.1  Introduction 

In  this  chapter  an  optimization  procedure  is  used  together  with  the 
analysis  program  described  in  the  previous  chapter  to  treat  certain  optimi¬ 
zation  and  design  problems  involving  arrays  of  parallel  thin-wire  antennas. 

As  mentioned  earlier  these  problems  include  design  of  unequally  spaced  arrays 
to  provide  equal-sidelobe  patterns,  reduction  of  pattern  sidelobes  through 
adjustments  of  wire  .lengths,  and  selection  of  wire  lengths  and  spacings  to 
optimize  array  directivity.  The  optimization  procedure  used  is  due  to 
Rosenbrock  [12]  and  represents  a  modified  version  of  the  method  of  steepest 
descent.  This  method  does  not  require  calculation  of  any  derivatives  and 
is  particularly  well-suited  to  automatic  calculation  by  digital  computer. 

The  method  has  been  used  by  other  antenna  engineers  in  treating  challenging 
optimization  problems  of  practical  interest. 
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of  Unequally  Spaced  Arrays 


The  analysis  program  presented  in  the  previous  chapter  has  been  used 
with  the  optimization  method  of  Rosenbrock  [12]  to  develop  a  design  technique 
for  linear  arrays  of  unequally  spaced,  parallel-wire  antennas.  As  before  the 
wires  are  assumed  to  be  centerfed  and  all  of  the  same  radius.  In  this  case 


it  is  also  assumed  that  each  wire  is  excited  with  a  real  unit  voltage,  and 
that  the  interelement  spacings  are  symmetrical  abov .  the  array  center. 


During  the  course  of  this  investigation  several  error  criteria  were 
tried  in  association  with  the  optimization  procedure  used.  Most  of  these 
resulted  in  relatively  poor  or  slow  convergence  to  acceptable  solutions,  and 
quite  often  an  acceptable  solution  was  never  reached  within  the  reasonable 
time  constraints  placed  on  the  computer.  Of  course,  the  error  criterion  is 
the  quantity  to  be  minimized  as  the  optimization  procedure  progresses  and  the 
most  successful  of  those  tried  is  given  by 


e  -  I  |T(1) 1 2  -  |T(i)|2  (8) 

i«2  1  1 

where  T(i)  is  the  peak  of  the  ith  sidelobe  with  T(l)  being  the  peak  of  the 
sidelobe  closest  to  the  main  beam.  (Note  that  the  resulting  broadside  pat- 
term  is  symmetrical  about  the  main  beam  since  the  excitations  are  uniform 
and  the  spacings  are  symmetrical  about  the  array  center.)  The  limit  p  is 
simply  the  number  of  sidelobes  to  be  controlled  on  either  side  of  the  main 
beam.  The  starting  point  for  the  iterative  procedure  is  the  pattern  of  an 
equally  spaced  uniform  array  with  a  given  number  of  elements.  The  result 
is  an  equally  excited,  unequally  spaced  array  with  nearly  equal  sidelobes. 
Of  course  the  design  parameters  are  the  interelement  spacings,  although  the 
spacing  with  respect  to  the  array  center  of  the  outermost  symmetrical  pair 
of  elements  is  kept  fixed  at  its  initial  value  so  that  the  bearawidth  will 
remain  essentially  unchanged  from  its  value  at  the  start. 

As  an  example  consider  an  8-element  linear  array  of  centerfed,  half¬ 
wave,  parallel  wires  that  are  half-wavelength  spaced  and  all  of  radius 
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0.007022A.  With  equal  in-phase  excitations  the  program  of  the  previous 
chapter  can  be  used  easily  to  obtain  the-  current  distributions,  input 
impedances,  and  the  H-plane  pattern  of  Fig.  7.  The  program  can  then  be 
used  in  conjunction  with  the  optimization  procedure  and  error  criterion 
given  by  (8)  to  obtain  a  new  pattern  with  reduced  and  nearly  equal  side- 
lobes  using  only  the  interelement  spacings  as  design  parameters.  The 
positions  of  the  elements  of  the  outermost  pair  are  kept  fixed  to  provide 
roughly  the  beamwidth  of  the  uniform  array  and  the  remaining  three  inter¬ 
element  spacings  are  adjusted  to  achieve  nearly  equal  sidelobes  at  a 
reduced  level.  The  resulting  pattern  is  also  shown  in  Fig.  7  along  with  the 
initial  and  final  spacing  configurations.  Note  that  because  of  symmetry  the 
spacings  are  shown  for  only  half  the  array  and  the  patterns  are  shown  for 
only  one  quadrant. 

Figures  8-10  show  results  of  applying  the  same  procedures  to  linear 
arrays  of  9,  12,  and  18  parallel,  half-wave  wires  respectively.  It  is  noted 
that  there  is  no  substantial  beamwidth  deterioration  even  though  the  sidelobes 
have  been  reduced  in  overall  level.  These  computations  were  all  performed 
using  an  IBM  370/155  computer  with  total  times  of  25,  36,  70,  and  185  seconds 
required  for  obtaining  the  results  shown  in  Figs.  7-10,  respectively. 


Finally,  it  should  be  noted  that  a  different  error  criterion  was  used 
earlier  in  conjunction  with  an  alternative  analysis  procedure  and  an  opti¬ 
mization  method  due  to  Fletcher  and  Powell  [13]  to  treat  the  same  problem 
discussed  above  [14],  However,  the  rate  of  convergence  proved  to  be  a 
serious  problem  in  that  case  and  the  present  choices  appear  to  be  much  more 
satisfactory  for  this  particular  problem. 


3-3  Use  of  Unequal  Wires  Lengths 


A  second  problem  treated  successfully  using  the  optimization  method  of 
the  previous  section  involves  use  of  array  wire  lengths  rather  than  the  inter¬ 
element  spacings  as  design  parameters.  The  objective  is  to  reduce  the  side¬ 
lobes  of  an  equally  excited,  equally  spaced  array  of  (initially)  equal-length 
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Fig.  7  -  Field  patterns  of  an  8-element  uniform  array  along 
with  the  pattern  of  a  similar  array  with  unequal 
spacings  designed  for  equal  sidelobes. 
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Fig.  8  -  Field  pattern  of  a  9-element  uniform  array  along  with  the 

pattern  of  a  similar  array  with  unequal  spacings  designed  for 
equal  sidelobes. 
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Fig.  9  -  Field  pattern  of  a  12-element  uniform  array  along  with  the 
pattern  of  a  similar  array  with  unequal  spacings  designed 
for  equal  sidelobes. 


wires  by  adjusting  the  wire  lengths.  As  before,  it  is  assumed  here  that  the 
resulting  array  remains  symmetrical  about  its  center  and  also  with  respect 
to  the  plane  containing  the  feed  points.  The  wires  are  again  assumed  to  be 
centerfed  and  all  of  the  same  radius. 

In  this  case  the  most  success  was  achieved  using  an  error  criterion 
given  by 

e,  -  l  1 T (i) | 2  -  |T  I  2  for  | T (i) |  >  |T  | 
x  i-1 

«  0  for  | T (i)  |  <_  |t  |  (9) 

where  in  this  case  T(i)  is  the  value  of  the  field  pattern  in  the  ith 
specified  direction  and  n  is  the  number  of  observations  calculated 
(e.g.  one  every  two  or  three  degrees)  over  the  sidelobe  region  of  the 
pattern.  The  quantity  j Tq | 2  is  simply  the  square  of  the  desired  sidelobe 
level.  The  starting  point  for  the  iterative  procedure  is  the  pattern  of 
an  equally  spaced  uniform  array  with  a  given  number  of  equal-length  wires. 
The  result  is  an  equally  excited,  equally  spaced  array  of  wires  of  dif¬ 
ferent  lengths  with  sidelobes  near  the  desired  level. 

As  an  example  consider  a  12-element  array  of  centerfed,  parallel 
wires  that  are  half-wavelength  spaced,  uniformly  excited,  and  initially 
all  one  half-wavelength  long.  The  wires  are  all  of  radius  0.007022X,  and 
the  initial  H-plane  pattern  is  shown  in  Fig.  9.  The  analysis  program  of 
the  previous  chapter  can  be  used  together  with  the  optimization  procedure 
and  error  criterion  given  by  (9)  to  obtain  the  new  pattern  of  Fig.  11  by 
adjusting  only  the  wire  lengths.  The  desired  level  in  this  case  was  +0.1 
corresponding  to  | T^  j  2  *  0.01.  This  corresponds  to  a  desired  sidelobe  level 
of  -20  dB.  Figure  11  also  includes  the  final  wire  lengths  which  are  assumed 
symmetric  about  the  array  center. 

Results  for  thi3  example  are  quite  close  to  those  desired  as  indicated 
by  Fig.  11.  However,  the  procedure  is  time-consuming  with  convergence  much 


slower  than  in  the  examples  of  the  previous  section.  A  total  of  about  nine 
minutes  was  required  to  obtain  the  data  corresponding  to  Fig.  11  using  an 
IBM  370/155  computer. 

3.4  Adjustments  of  Both  Lengths  and  Spacings 


An  obvious  extension  of  the  work  for  the  last  two  sections  involve 
adjustments  of  both  wire  lengths  and  interelement  spacings  simultaneously 
to  achieve  reduced  sidelobes.  Figure  12  shows  the  final  result  of  mini¬ 
mizing  the  error  criterion  given  by  (9)  for  an  8-element  array  of  parallel 
wires.  The  starting  point  for  the  iterative  procedure  was  the  initial 
pattern  of  Fig.  7  corresponding  to  a  uniform  8-element  array  of  half-wave 
wires  that  are  half-wavelength  spaced,  and  the  final  result  is  that  of 
Fig.  12  while  Fig.  13  indicates  the  corresponding  array  configuration. 
Symmetry  about  the  array  center  and  about  the  plane  containing  the  feed 
points  is  again  assumed  and  the  wire  radii  are  all  0.007022A.  The  result 
is  quite  good  although  convergence  was  considerably  slower  than  in  the 
case  involving  just  the  spacings.  The  results  of  Fig.  12  required  about 
11  minutes  using  an  IBM  370/155  computer. 

As  a  final  illustration  of  the  use  of  the  analysis  program  described 

in  Chapter  2  together  with  Rosenbrock’s  optimization  procedure  consider  the 

design  of  a  typical  Yagi  antenna.  Here,  lengths  and  spacings  are  to  be 

determined  for  an  array  of  parallel  wires  (with  symmetry  assumed  only  about 

the  center-points  of  the  array  wires)  such  as  to  maximize  the  endfire  gain 

G  defined  as 
o 

m  4tt  (radiation  intensity  for  the  endfire  direction) 
o  power  input  to  the  array  '  ' 

Formulas  for  computing  the  gain  using  matrices  available  from  the  analysis 
program  are  well-known  [1,2],  These  can  be  used  with  the  iterative  opti¬ 
mization  procedure  to  optimize  a  Yagi  antenna  with  respect  to  its  wire 
lengths  and  spacings.  The  starting  point  can  be  an  initial  state  correspond¬ 
ing  to  some  classical  or  other  improved  design  procedure.  Figure  14  shows 
results  after  three  iterations  for  Yagi  antennas  numbering  3-6  elements  where 
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Fig.  12  -  Pattern  of  an  8-element  array  with  sidelobe  level 
minimized  with  respect  to  both  lengths  and  spacings 
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in  each  case  the  initial  state  (already  the  result  of  some  gain-maximizing 
design  procedure)  was  at  least  slightly  improved  upon.  For  the  3-element 
case  the  initial  state  was  the  result  of  a  classical  design  procedure  [15] 
and  in  the  other  cases  the  initial  status  were  due  to  another  analysis  and 
optimization  program  based  on  matrix  methods  [16].  The  time  required  to 
optimize  the  6-wire  Yagi  antenna  was  2  minutes  on  the  IBM  370/155  computer. 
Improvements  in  gain  were 

a)  from  5.9  to  9.1  for  3-element  case,  starting  from  an  antenna 
designed  by  a  standard  technique  [15] 

b)  from  14.2  to  14.5  for  4-element  case,  starting  from  an  antenna 
optimized  by  matrix  methods  [16] 

c)  from  18.9  to  19.4  for  5-element  case,  starting  from  an  antenna 
optimized  by  matrix  methods  [16] 

d)  from  24.4  to  24.9  for  6-element  case,  starting  from  an  antenna 
optimized  by  matrix  methods  [16]. 

The  examples  presented  in  this  chapter  are  by  no  means  exhaustive 
and  serve  simply  to  illustrate  possible  uses  of  the  analysis  program  de¬ 
scribed  in  Chapter  2  and  presented  in  Appendix  B.  The  examples  do  point 
out,  however,  that  this  particular  analysis  program  can  be  used  in  solving 
quite  difficult  problems  without  requiring  unreasonable  computer  time  as 
is  often  the  case. 


4.  A  SPECIAL  PROGRAM  FOR  HANDLING  CERTAIN  WIRE  GROUND  SYSTEMS 

4.1  Introduction 

'  < 

This  chapter  describes  a  special  user-oriented  computer  program  that 
was  designed  specifically  to  handle  certain  problems  involving  vertical  wire 
antennas  over  radial  ground  systems.  The  program,  presented  in  Appendix  G 
of  this  report,  is  capable  of  treating  fairly  general  radiation  problems  for 
wire  configurations  such  as  depicted  in  Fig.  16.  The  program  is  a  speciali¬ 
zation  of  the  general  analysis  program  presented  in  an  earlier  Scientific 


Report  [4]  that  is  capable  of  handling  radiation  problems  involving  arbitrary 
configurations  of  thin  wires  including  junctions.  By  incorporating  available 
symmetries  the  program  presented  here  offers  considerable  advantage  over  the 
arbitrary  wire  program  with  regard  to  both  speed  of  computation  and  simplicity. 

The  typical  problem  geometry  is  assumed  to  consist  of  a  z-directed 
center  wire  plus  a  number  NB  of  branches  as  in  Fig.  15.  The  branches  are 
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assumed  to  be  all  the  same  length  and  all  at  the  same  angle  with  the  center 
wire.  The  branches  are  uniformly  distributed  about  the  center  wire,  so  that 
they  are  situated  at  regular  intervals  of  2tt/NB  radians. 

As  in  the  previous  program  all  wires  are  assumed  here  to  be  perfect 
conductors.  Each  wire  length  L  and  radius  a  are  such  that  L/a  »  1  and 
a  «  A.  Within  Harrington’s  general  procedure  the  so-called  triangle  func¬ 
tions  are  used  for  subsectional  current  expansion  functions  resulting  in  a 
piecewise  linear  approximation  to  the  current  of  each  wire.  Triangle  func¬ 
tions  are  also  used  for  testing  or  weighting  functions  resulting  in  a 
Galerkin  solution  co  the  analysis  problem  which,  as  mentioned  earlier,  is 
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characterized  by  relatively  fast  convergence.  The  wires  are  assumed  to  be 
unloaded  and  the  possible  locations  of  applied  excitations  are  restricted 
to  points  along  the  z-directed  wire  (see  Fig.  15)  which  correspond  to  peaks 
of  triangle  current  expansion  functions.  Most  of  the  necessary  details  of 
the  program  were  presented  earlier  with  the  arbitrary  wire  programs  [A].  The 
only  details  included  here  are  those  that  differ  significantly  from  the 
previous  work. 

4.2  Data  Input 

All  input  data  are  provided  through  the  main  program.  The  first  data 
statement  reads  in  quantities  NB,  NEB,  NEC,  NR,  BK,  BLENTH,  CLENTH,  and  ANGLE 
where 

NB  »  the  number  of  branches  in  the  wire  configuration  not  including  the 
center  wire.  For  example,  in  Fig.  15a  NB  *  2  while  in  Fig.  15b 
NB  *  4. 

NEB  =  the  number  of  triangle  expansion  functions  chosen  for  each  branch. 

(All  branches  have  the  same  length.  The  center  wire  is  not  included 
among  the  branches.) 

NEC  *  the  number  of  triangle  expansion  functions  chosen  for  the  center  wire. 
(It  should  be  noted  with  regard  to  both  NEC  and  NEB  that  for  accurate 
descriptions  of  current  distributions  it  is  advisable  to  use  at  least 
12  triangle  current  expansion  functions  per  wavelength  of  wire.  On 
the  other  hand  a  smaller  number  may  be  used  if  only  far-field  patterns 
are  of  interest.) 

BK  =  the  factor  2ir  divided  by  the  wavelength  in  meters. 

BLENTH  =  the  length  in  meters  of  each  branch.  (As  mentioned  earlier, 
all  branches  have  the  same  length.) 

CLENTH  *  the  length  in  meters  of  the  center  wire. 

ANGLE  *  the  angle  between  the  branches  and  the  positive  z-axis.  (All 
branches  are  assumed  to  be  at  the  same  angle  with  the  center 
wire.)  This  angle  is  expressed  in  degrees  so  that  in  Fig.  15a-c 
ANGLE  =  90  and  for  Fig.  15d  ANGLE  -  180. 


The  second  data  statement  provides  the  wire  radius  in  meters.  (For  the 

example  included  with  the  program  cf  Appendix  C  the  wire  radius  is  given 
-L 

as  10  ‘  meters.)  The  third  reads  in  a  number  NF  which  equals  the  number 
of  independent  feed  voltages  to  be  applied  to  the  center  wire.  As  mentioned 
earlier  these  feed  voltages  must  be  applied  at  wire  positions  corresponding 
to  peaks  of  triangle  current  expansion  functions  so  that  NF  must  not  exceed 
the  number  of  triangles  assigned  to  the  center  wire.  For  radiation  problems 
it  is  assumed  NF  >_  1,  and  DO  LOOP  44  executes  a  total  of  NF  times  with  an 
integer  J1  and  a  complex  number  V  read  in  using  a  single  additional  data 
card  with  each  execution.  J1  is  the  number  of  the  particular  triangle 
whose  peak  marks  the  desired  position  of  the  first  excitation  voltage  and 
V  is  the  desired  value  of  that  excitation  in  volts.  Thus,  following  the 
data  statement  providing  NF,  the  next  NF  data  cards  read  in  feed  positions 
(corresponding  to  the  numbers  of  specific  triangles)  and  excitation  voltages 
for  each  of  the  feed  points  along  the  center  wire. 

Far-field  patterns  are  calculated  and  printed  out  using  DO  LOOP  17. 
Here,  of  course,  it  is  necessary  to  provide  the  polar  and  azimuth  angles 
of  observation  desired.  In  the  example  of  the  program  included  in  Appendix  C 
the  instructions 

DO  17  IPH  -1,2 
PHI  =  (IPH-1)*  90 
DO  17  ITH  =  1,181,5 
THE  -  ITH-1 

cause  field  patterns  to  be  computed  in  five  degree  steps  of  the  polar 
angle  9  for  two  different  choices  of  azimuth  angle  <j>  =  0° ,  90° . 


Program  Description 


Details  of  the  problem  geometry  are  calculated,  from  the  input  data 
using  subroutine  GOMTRY  and  the  generalized  impedance  matrix  [Z]  is  com¬ 
puted  using  subroutine  ZBRCH.  [Z]  is  a  square  matrix  of  dimension  equal  to 
the  number  of  independent  expansion  functions  used,  and  its  elements  are 
functions  only  of  the  problem  geometry.  This  matrix  is  inverted  using  sub- 


routine  LINEQ  and  the  current  matrix  [ I ]  is  calculated  by  using  this  inverted 
matrix  in  (1)  and  DO  LOOP  9.  Finally,  the  field  patterns  specified  are  com¬ 
puted  and  printed  out  by  way  of  DO  LOOP  17  as  mentioned  earlier.  In  passing 
it  should  be  pointed  out  that  in  the  program  presented  in  Appendix  C  the 
matrix  [Z]  1  occupies  the  same  storage  locations  previously  held  by  [Z]. 
Hence,  the  term  [Z]  denotes  the  generalized  impedance  matrix  before  inver¬ 
sion,  and  the  generalized  admittance  matrix  [Z]  ^  after  inversion.  It  should 

also  be  noted  that  [Z]  is  manipulated  in  this  program  as  a  column  matrix  of 
2 

N  complex  numbers  where  N  is  the  number  of  independent  expansion  functions 
used. 

4.3  Examples 

To  illustrate  use  of  this  special  purpose  program  consider  the  wire 
configurations  shown  in  Fig.  15.  For  the  antenna  of  Fig.  15a  suppose  the 
length  of  the  center  wire  and  the  lengths  of  the  branches  are  all  A/4. 


Fig.  16  -  Junction  treated  as  overlapping  open-ended  wires. 

Suppose  further  that  the  wavelength  is  given  to  be  one  meter.  Thus,  BK  *  2ir 
CLENTH  *  0.25  and  BLENTH  =  0.25.  In  accordance  with  previous  discussions  of 
the  treatment  of  wire  junctions  [5]  the  antennas  of  Fig.  15  can  be  treated  as 
systems  of  overlapping  open-ended  wires  as  shown  in  Fig.  16,  corresponding  to 
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the  antenna  of  Fig.  15a.  As  indicated  the  antenna  can  be  analyzed  con¬ 
veniently  using  six  triangle  current  expansion  functions  for  the  center 
wire  and  seven  for  each  branch,  so  that  NEB  *  7  and  NEC  *  6.  Obviously, 
the  branches  are  normal  to  the  center  wire  in  this  case  so  that  ANGLE  *  90. 

The  radius  of  each  wire  is  0.007X  so  for  the  second  data  statement  RAD  *  0.007. 
Finally,  the  antenna  is  assumed  to  be  base  fed  with  a  real  unit  excitation 
voltage.  Hence,  NF  *  1  and  since  the  base  corresponds  to  the  peak  of  the 
seventh  triangle  the  fourth  data  statement  provides  J1  *  7  (locating  the  feed 
point)  and  an  excitation  given  by  V  ■  1.0  +j0.0  volts. 

If  the  wire  lengths  and  radii  are  the  same  as  in  Fig.  15a  then  the 
antennas  of  Figs.  15b-d  can  be  handled  with  very  simple  changes  in  the  input 
data.  For  example,  if  the  antenna  of  Fig.  15b  has  a  center  wire  and  four 
branches,  all  of  length  A/4,  and  if  all  remaining  data  are  the  same  as  before 
the  only  change  required  in  the  input  is  NB  **  4.  The  vertical  wire  is  always 
taken  to  be  z-directed  and  the  projection  of  the  first  branch  on  the  x-y  plane 
is  always  assumed  to  be  x-directed,  with  the  projections  of  the  remaining 
branches  distributed  at  equal-angle  intervals.  Thus,  in  Fig.  15b  the  first 
branch  assumes  the  positive  x-direction,  the  second  the  negative  y-direction, 
and  so  on.  Next,  for  the  L-shaped  antenna  of  Fig.  15c,  it  is  obvious  that  the 
only  change  in  input  data  required  is  NB  *  1  if  all  other  characteristics  are 
the  same  as  in  the  first  example.  Of  course,  in  this  case  the  only  branch  is 
+x~directed.  Finally,  the  straight  A/2-wire  of  Fig.  15d  can  be  treated  as  a 
A/4  z-directed  wire  with  a  single  branch  of  length  A/4  at  an  angle  of  180°  with 
the  +z-axis.  Here,  input  data  differing  from  the  specifications  of  Fig.  15a 
are  NB  =  t  and  ANGLE  =  180. 

Current  distributions  for  the  antennas  of  Fig.  15  are  plotted  in 
Fig.  17.  Input-output  data,  current  distribution,  and  sample  field  patterns 
for  an  additional  example  are  included  with  the  program  in  Appendix  C. 


3 


Fig .  17a  -  Current  for  Fig.  15a. 


Fig.  17c  -  Current  for  Fig.  15c. 
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5.  CONCLUSION 

Two  user-oriented  computer  programs  have  been  presented  and  described 
for  handling  radiation  problems  involving  certain  special  configurations  of 
thin-wire  antennas.  Both  of  these  stem  from  the  method  of  moments  with  the 
first  employing  piecewise  sinusoidal  current  expansion  functions  and  the 
second  a  piecewise  linear  current  approximation.  The  first  treats  radiation 
by  arrays  of  thin,  parallel,  centerfed  wires  where  the  feed  points  all  lie 
in  the  same  plane.  The  wires  are  assumed  to  be  lossless  and  all  of  the  same 
radius.  To  illustrate  use  of  the  program  several  array  design  problems  were 
treated.  These  include  design  of  unequally  spaced  arrays  to  provide  equal- 
sidelobe  patterns,  selection  of  wire  lengths  and  spacings  to  optimize  direc¬ 
tivity,  and  reduction  of  pattern  sidelobes  through  adjustments  of  wire 
lengths . 

The  second  program  presented  here  treats  efficiently  certain  special 
situations  involving  vertical  wire  antennas  over  radial  wire  systems.  Here 
again  the  wires  are  assumed  lossless  and  all  of  the  same  radius,  and  feed 
points  are  restricted  to  the  vertical  wire.  Several  examples  were  included 
to  illustrate  use  of  this  s*:und  program  as  well. 
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APPENDIX  A 

EVALUATION  OF  THE  ELEMENTS  OF  THE  GENERALIZED  IMPEDANCE  MATRIX  [Z] 


For  an  array  of  straight,  z-directed  wires  the  elements  of  [Z]  are 
given  by  (6)  and  can  be  evaluated  as  follows: 


E  •  dz 
m 


Z  2  " 

n  sin  k(z-z  ,)  t  n+1  sin  k(z  ,,-z)  ,_n 

_ _ _ sdL  +  _ L-StL-l  x _ 152 _ 

sin(kAz  )  I  sin (kAz  )  sin(kAz  ) 


-jkR1  -jkR2  -jkR3 

- - 2  cos  (kAz  )  - - +  - - 

ra  *^2 


where 


Ri  *  /a*  +  (z-z  .)' 
1  m-1 


R.  -  /a*  +  (z— z)  ^ 
z  m 

R3  ’  ^  +  ('-VlV 


Az  *  z  -z  ,  *  z  ,,-z 
n  n  n-1  n+1  n 


Az  =  z  -z  ,  *  z  ,,-z 
m  m  m-1  m+1  m 


The  computed  elements  are:  Z  «*  R  +  IX  where, 
r  nm  nm  J  nm  ’ 


P  m  -  -  -  ,  -  -  -■  _ _  x 

nm  sin(kAz  )  sin(kAz  ) 
m  n 


[cos  kCz^  -  (CKvq)  +  Ci(ug)  -  CKi^)  -  CKv^} 


+  sin  lc(zia_1  -  zn_^)  (Si(Vg)  -  Si(Ug)  +  SiCu^)  -  Si(v^)} 


+  cos  -  z^_^)  {Ci(v^)  +  Ci(u^)  -  Ci(Ug)  -  Ci(v^)} 


+  sin  ~  z  j.)  {Si(v^)  -  Si(u^)  +  Si(u5>  -  Si(v5)} 


-2  cos(kAz  )  cos  k(z  -  z  ,  (Ci(v0)  +  Ci(u0)  -  Ci(u_)  -  Ci(v  )} 
m  m  n-1  2  2  3  3 


-2  cos(kAz Jcos  k(z  -  z  .)  Si(v0)  -  Si(u„)  +  Si(u_)  -  Si(v.)} 
th  m  n-1  l  i  j  j 


+  cos  k(zm_1  -  zn+1)  {Ci(Vg)  -  Ci (v^)  +  Ci(u^  -  Ci(u1>} 


+  sin  k(z^_^  -  z^+1)  {Si(Vg)  -  Si(ug)  +  Sifa^)  -  Si(v^)} 


+  cos  kU^  -  zn+1)  {Ci(vQ)  -  Ci{v,.)  -  Ci(u-)  +  Ci(uR)} 


+  sin  k(znrt.1  _  zn+1)  tSi(Vg)  -  Si(Ug)  +  Si(u5>  -  Si(v5>} 


-2  cos(kAz  )  cos  k(z  -  z  .)  {Ci(v7)  -  Ci(v,)  -  Ci(u_)  +  Ci(u7)} 
m  m  n+l  /  j  j  / 


-2  cos(kAz  )  sin  k(z  -  z  .)  {-  Si(u7)  +  Si(v_) 
m  in  n+l  /  / 


+  Si(u3)  -  Si(v3)}] 


'«MwoUru  KUMMm  KM'*  •*** 


A<3*(  !-l)+J)*ZN+7* 

*>  7y*ZM+nzM 

4  7N=ZN*nZN 

h  00  3  1*1 • 9 

IF  (All)  «LT •  0.0)  GO  TO  21 
'MI)*S0PT(RR4-Ain*Mn)+M  n 
VII )=RR/U< I ) 

GC  TO  22 

?1  V ( I )  =  SOPT  CRR+AC I)*A(I))-A(I) 
till  )*RR/VU  ) 

2?  CALL  SICKSI  ,CI,U(  !)> 

six  n=si 
ax  n-ci 

call  tchsi.ci  fv<  m 

SV!  I )  =  S  I 

cv ( n  =c  i 

sp{  i  »=su(  n+svi  i ) 
sm  n=su<  i  »-sv<  i ) 

C P( I  )*CU( i )+cv( I  ) 

CN(  I  )=C'I(  I  )-CV(  I  ) 

’  Cf'NT  INIIF 

cm*cosu<  n ) 

r.(2)=cns(A(2) » 

C(3)  =  C0S(ACM  | 

C(7)  =  cns( A( 7)  ) 

C(8)*cns( A( R) ) 

C(9)=CnS(A(Q) ) 

SI  I )  =  SIN(A( 1) ) 

S I  2  )*S INI A( 2 )  ) 

sm=sTNU(  ■*) ) 

SI7)*SIN(A( 7) ) 

Sm-S1N(A(P)  ) 

S(9)=SfN(A(Q) ) 

Pl=C(  1>*(CP(1)-CP(4M-S(  l)*(SN<4l-SN(U  )4-C(  3)*f  CP(3»-CD(6)  ) -S  t  ?)* 
f  (SN{<S  )-SN(*|  t+C(  7»*(CP(7)-CP(4))-S(  7)  *1  SN  ( 4  l-SN(  7  U  +C.  ( 9 )  *  I C  0  ( 9  ) 
$-CP(6)  )-S«9)*(  SN(  f  I-SM9)  ) -COM  Cl?)  MCP 1 2  l-CP  I  •>  I) -S  l  2 )  *  I SN  (  5  ) 

*  -SN(2tH-C(ft)*(CPf  Bl-CPfSI  l-SIBIM  SNfSI-SN(fi)  )  ) 

AG=C(  1)*(  SPIM-SPI  l)  )-S(  1)*ICN(4)-CN(  l)  )+C(  3I*(SPIM-SP(3)  )-S(  3)* 

♦  (CN(6l-CNm  )+CI7i*(SP(4)-SP<  7)  )-S( 71 *<CN(4 )-CN(7 It +C ( 9  I *(  SP(*  )- 
*SP(<?)  )-$(9)*(CN(6)-CN(9))-CC*(C(2)*(SP(5)-SP(?)  )-.?(?)*  (CNC  5  )-CM 
S?))*C(*)*(SPf  5J-SPIR)  )-S(R)*<CNf5l-CN(*))  ) 

7MN=CMPLX(RL»  AG)*15./(SIN(0ZN»*SIN(0ZM)  ) -*- 7«N 
I F f N.FQ.2  .OP.  N.F0.3)  RFTURM 
GO  TO  10 
FNO 


ono  o  o  o  ooo  r>  o  o 
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APPENDIX  B 

PROGRAM  FOR  LINEAR  ARRAYS  OF  PARALLEL  WIRES 

This  program  is  suitable  for  analyzing  radiation  by  linear 
arrays  of  parallel  thin-wire  antennas.  The  wires  can  be  of  unequal 
lengths  although  all  wire  radii  are  assumed  to  be  the  same.  Sym¬ 
metry  is  assumed  only  about  the  plane  containing  the  midpoints  of 
the  wires.  The  wires  are  assumed  to  be  lossless  with  no  externally 
applied  loading.  The  sample  output  corresponds  to  analysis  of  an 
8-element  array  of  half-wave  wires  that  are  half-wavelength  spaced 
and  all  of  radius  0.01A.  Five  expansion  functions  are  used  for 
approximating  the  current  of  each  wire.  (Because  of  symmetry  only 
three  of  these  are  independent.)  Each  wire  is  fed  with  a  real  unit 
excitation  and  the  principal  H-plane  pattern  is  computed  in  steps 
of  three  degrees. 

The  program  is  described  in  Chapter  2.  The  main  program 
is  presented  first  followed  by  the  subroutines. 


***  MAIN  PROGRAM  *** 

DIMENSION  V C  20 1 ♦ VV( 10 ) 

COMMON  AK ,N,NL,N7,NV*VV»MQ 
COMPLEX  CMPLX,CFXP,VV 
TP=6.2«31R53 

A«  =  RADIUS  *  <  2*P I  ) 

AK=0.01*TP 

N  =  TOTAL  no.  of  elements  in  the  linear  array 


N  =  P 

NV  =  TOTAL  NO.  UE  VAR  I  ARIFS  l  WIRE  LENGTHS  f.  I  NOE  PENDENT  ARRAY  SPACING*;  > 


NV=  15 
C 

C  FIRST  FFW  ELEMENTS  PE  YU)  DENOTF  THE  LENGTHS  OF  THE  ELEMENTS  AND 
C 

C  THE  REST  DENOTE  THE  DISTANCES  OF  THE  ELEMENTS  FROM  THE  FIRST,  WHICH 
C 

C  IS  CONSIDERED  AS  THE  ORIGIN  OF  THE  ARRAY. 


oooooo  o  o  o 


c 

C  ***  NOTE  ***  AIL  DIMENSIONS  OF  Y( I )  ABE  METFRS/WAVELENGTH  *  (7*PII 

C 

C 

Y ( 1 )=0.2?*TP 
Y(2)=0.25*TP 
Y(»)=0.2«*TP 
Y(4)=0.2B*TP 
Y(5)=0.25*TP 
Y(M  =  0„25*TP 
Y(7)»0.2‘5*TP 
Y(fl)=0.25*TP 
Y{9 )=0.50*TP 
v(  1D)=1.0*TP 
Y ( 1 1 ) s 1 • 5*TP 
Y ( 1 2 )  =  2  »0*T  P 
Y (  n)=2.5*TP 
Y( 14)=3.0*TP 
Y(  1M=3.*?*TP 
C 

C  SPFCIFIFS  THE  FEED  VOLTAGES 

C 

OP  3  1=1 tN 

VV(  n=CMPLX(1.0,0.0) 

■»  CONTINUE 
C 

C  PHl'CIPAl  H  -  PLANE  PATTERN  IS  COMPUTFP  IN  STFPS  OF  MQ  DEG»EFS 


MO=3 


*****  FND  OF  DATA  ***** 

N7=N+1 

Ml=N-l 

CALL  ASCTFO(Y) 

STOP 

END 


StlPROiJTTNE  ASCTFDI Yl 

Cr«PLFX  Z»V,SS,DETEPM,CEXP,CMPLXtTERMf A7,ZMN»YJ,7M1 ,yV 
COMMON  AK,N,NL»NZfNVtWtMO 

DIMENSION  N(15)fYZ(10)»n7(l5ltZ(50*S0I#V(S0»l)»SS(50»l)»CSnSI 
DIMENSION  TE(200)tPHl200),Y{20»  tVV( 10) 

B 10=0.0 

Min=o 

DO  2  I  =  1 » N 


*****  NOTF  ***** 

FOR  OPTIMIZATION  PROBLEMS  MM=  INT I  YID-O  .  ST1 )  *1 

FOR  ANALYSIS  PROBLEMS  mm«I NT(Y( I I-0.5T1)*3 


m( 'I  ♦!»*»(  I ) 

H7(  I  »*V(  !  I /TLOATl **M| 

MK;sM  f  v+  l ) 

■V-  6<-  1“  1 « MK 

VI  T  .1  x(-0.0#0.0) 

on  4  j*.it* 

M??  = MI.D  +  l 

yc«*2?  •  i-i*wf  .i » 
n  cm  in.i  ip 
¥7(11=0.0 
or  s  J*  i  ♦  ni 

Y7I.Ui»  =  Y(M  +  J) 
on  fit-  i  =  I  *  n 
7N=n.O 
«ll*M(Llfl 
Ml  2* M  (  I  *\  ) 
nn  SI  ! =ml 1 . ml 2 
7M*0.0 
nn  S? 

H  I  tJl=7Mll  7v,7MtnZ(U,AKI 
7.^=7M  +  f:7(L) 

2f'S/.\.  +  P7(L) 

rC-M  !.•-'(  if 
no  74  k  =  l.Ni 
K1=K+1 
7i'r-r.7  Ik  ) 

MK  l  (  K  )  + 1 
MK?*Y(K  +  1-| 

OP  71  !sMKttM«2 

nn  7^ 

z Yr-H7 (  I  ) 

Ml  l«H(UM 
Ml ?  *M ( l 4 1 J 

nn  f>  j*ml i tMi? 

SD=Y7IU-YZ(K) 

71  I  tJ)«7Mfc(7M»?M»07lLI»r)7(K|.SD) 
7Mr7M«-n7(L) 

rnNTTwuF 

7M=7M*07(K) 

r.UNTlMUF 

nr  A4  K  =  ?,N 

K1  =  K-1 

ZK>*-P7(M 
MN1=M(K  Ul 
fc^2*M(K  +  H 

nn  fti  i=vMifMN2 
nn  ft?  I «l,Kl 
7M=-n zn  ) 

Ml l*M(ll+l 

Ml?«M(l+l \ 

on  ft2 

SD*YZ(K  )-Y7(l  ) 
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7II*J|*7MN(7M#?N,nZlL»#OZlKl,Sf)l 

62 

63  CONTINUE 

61  7N»2N+r»7<KI 

64  CONTINUE 

CALL  CSMTNIZ.MN, DETERM) 

CALL  MULfPYCMN#NN,l ,7*V.SS) 

LL*0 

no  11  LK*  1* T81 *MQ 
PC*0,6l744A3*Fi.0AtCLK-l ) 

«»c*rrisipc) 

A7*CMPLXC6. O.O.O) 

DO  44  J*1*N 
YJ*C>PLXI0.C*PC4Y7I  Jl ) 

TFP  M*f.  EXP!  YJ) 

CS<  J)*TAN(0.5*DZ!J)) 

MC*M(  J|  4-1 
M8=M( J+l) 

MAxMC.fl 

A7*AZ*SS( MC . 1 l*TERM*CS( J I 
inMA.GT.MR)  GO  TO  44 
on  45  I«MA,MR 

46  AZ* AZ.+  SSI  I .  1 ) *2.0*TERM*CSLJ ) 

44  CONTINUE 
Lt*n.+i 

TE ( LL  )*CABS I AZ I 

IFUElLL)  .GT.  RIG)  BIG*TE!LL) 

PHI LL )*ATAN2I AIMAGf AZI .REAL  I AZ) ) /O. 01  74533 
11  CONTINUE 

00  43  I»l,NV 
Y ( I  )*YII)/6. 2331853 
43  CONTINUE 

HP  I  TE  (  3*111 )  (Yd).  I«l,N) 

WRITFI3.112)  (YII).  I»N7,NV) 

112  FORMAT! •O*.*  SPACING’  ' .10F10.5/ //) 
ill  FORMAT! ' 0* »  *  LENGTH*  • ♦ 10F10. 5/// I 
WRITE!  3*892 ) 

892  FORMAT! »0». *  CURRENT  DISTRIBUTION  ON  THE  ELEMFNTS  •) 
00  49  K* 1  *  N 

MR l*M(K )4 l 
MR  p«M  ( K  ♦  T) 

WRITE! 3.895 ) 

WR I TE! 3. 888  )K 

883  FORMAT! *6*  *  *  CURRENT  ON  * . 1 5. 5X, * ELEMFNT* / ) 

00  50  NZ*MK1,MK2 
WRITE! 3.999 )  SS(N7,1) 

999  FORMAT*  •  S2F20.6) 

50  CONTINUE 
49  CONTINUE 

WRITF13.895) 

895  FORMAT  PO*  // ) 

WRITE! 3, 893) 

893  FORMAT! *0*. ’INPUT  IMPEDANCES  AT  THE  FEFO  POINTS'///) 
DO  51  J*1.N 

M22*M( j)4i 


•>  o  o  o  r»  o  r> 
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IFCCAPSIV(v?2,H).  EQ.  O.O)  GO  TO  51 
A7*V(M22,U/9SfM22,l) 

WRITF<3,R94|  J,AZ 

FORMAT |-«0 **  •' INPUT  IMPEDANCE  OF  1 » l 5  »5X» •  cLFMpNT  =«,?F20.6//» 

CONTINUE 

WPITF(3,B91 ) 

FORMAT  C • 1 • »  *  NntfMALTZED  F-  FIFLO  PATTFRN  «///) 

LMxO 

DO  Q5  t=i,u 
TFIL)«TFIl )/PlA 
wPITF(3,34?U  )  ,PH(L) 

FORMATC  *  ,  I5»?F40.M 

LFxLV+^O 

CONTINUE 

WE  ITf  ( 1.PQ0)  «IG 

FCRmAtMoN*  FORMAL  I  7  AT  I  ON  CONSTANT  ='.F15.5///I 

R  FT  URN 

FNP- 


FUNCTION  Z«N(7Vlf7NI,nZHtnZN,Rl 

CfMPUTt S  MUTUAL  IMPEDANCE  BETWEEN  TWO  ELEMENTS  oc  LENGTHS  07”  P7N 
c/ii-h  OF  RADIUS  »  HAVING  STARTING  POINTS  AT  7”1  AND  7NL 


***  NOTE  ***  All  DIMENSIONS  ARE  f*?*0!) 


C(”PLEX  7MN,CMPLX 

m  MENS  ION  A (9) ,U<9) tV<9) ,SUf9| ,SV(o) ,CUI9) ,CV<9),CC9) ,SJ9| ,SP(Q) , 
tSM{Q) fCN(Q) » CP  f  9 ) 

N  =  1 

IM7MI.F0.-0ZM)  N  =  2 
rC=2.C*CPS<P7M) 

kr=r*p 
7N=7N1 
DO  1  1=1,3 
ZM=7.M1 

nc  2  J  =  l»  3 

A  «  3* ( I— 1 1  ♦ J )  =  ZN-7M 

7N=ZM*P7M 

7N=7NfPZN 

7MN=CMP1  X(0. 0,0.0)  ><^| 

Gr  T°  6 
Z  N  =  7  N  l 

«*3 

PC'  4  1  =  1,3 
ZM=7.Ml 
OP  5  J*1 , 3 


^kAuL^ui  'a  L'tsim 


A 1 3*1  I-l)*J.)*ZNtZ»* 

•i  7y»z**h  z* 

4  ZN*ZN*OZN 

6  00  3  I * i»  9 

IF  (AMI  .LT.  0.0)  GO  T(1  21 
'i(  I  )*S0PT(RR+At  I  )*AU )  )  ♦Af  l  ) 

VI- 1  )*RR/UI  D 
GO  TO  22 

21  VI  I  )*SOPTCRR+AM  >*AM  U-AI  I  ) 

-tl't  I  )*RR/V(  If i 

2?  CALL  SICHSI.CI.UU  M 

SU( !)*$! 

Cl'M.)*CI 

CALL  TC 1 15  I «C! .Vi M ) 

svi  n*si 

C  Vi  I )  *C:I 

spi  n*sui  n+svi  i  > 
snuusumi-svim 
r.p(  i)*cuin>cvm 
CNm»cum-cvM  j 

*>  CONTINUE 

C 1 1 ) *COS(,AC  1) ) 

C(2)=COSIA(  2)  I 
C<  3)*C0S( A ( 3  V ) 

C(7)aC0S( A( 71 ) 
cm*cnsu(9) ) 

C(9)=CnS(A(P) ) 

sii)*sinu<  i) ) 

S(2)«SIN(A12)» 

U3)»SfN(Af ?) ) 

S(7)«SIN(A(7n 
S(8)*SIN(AIR) ) 

S(9).«;tN(A<0>  ) 

Rl*C(  l)*(CPm-CP(4))-SI  l)*(SN(4)-SNm  )+C(3)*ICP(3»-CP(6)  »-SI?)* 

f<SN(4)-SN(^|)+C{7l*(CP(7)-CP(4n-S(7)*ISM(4)-SN(7»)*C(9)*(CD<9J 
1— CP (6 ) )— S(9  )*( SNC6I“9N(9))-CC*ICI2)*CCP(2I-CPI5))-SI2»*(SN(5) 

*  -SNI2))+C|fl)*(CPI  8I-CPI5 SN(S)-SM(ai) ) 

AG*C(  l)*(SP(A )-SP(l))-S(ll*fCN(4)-CN(  l)  I4CI  3|*«$P(M-SP(3)  I-SI3)* 

♦  I  CM  (6  )“CN  (3  ))^CI7)*(SP(4)-SP(7))-SC7I*ICN(4)-CN(7))4-C(9I*(SPI4  )- 
HSP(9n-S(9)*{CNC6)-CN(9))-eC*IC(2)*CSP(5)-SP(?)l-SI2)*(CNI‘5)-rM 
12 ) l+Cl  *)*( SPI 5 )-SPI 4)  )-SI8)  *(CNI  51-CNCR ) )  ) 

7MN*CMPLX(RL,AG)*15./(S!N(0ZN)*SIMIDZM)  )«-Z*N 
IFIN.E0.2  .OP.  N.FQ.3)  RFTURN 
GO  TO  10 
FNO 


nonoonoa 


FUNCTION 


7*1(7.»«£N$0K9AM 


COMPUT&S  MUTUAL  IMPEDANCE  BETWEEN  TWO  ELEMENTS  OF  RADIUS  AK  AND 


LENf.TH.DM  HAVING  STARTING  POINTS  AT  ZM  AND  7N 
A**  NOTE  ***  ALL  DIMENSIONS  ARE  I*2*PM 


COMPLEX  CMPLX,ZM1 

cpk*cpsioki 

S  DK  *  l  S  *  /Tl  ,-CDK  A.f.DK  I 
n$Q«AK*AK 
p4=A.*cnk 
0D4  =  ?  C4*f  UK. 

n«i 

IFTZM.eo.C.CL  N=c 

A*ZN-ZM 

R*A-nk 

C=A+OK 

op -dm 
F*C*DK 
CA*COS( A  )■ 

S'A*  S  I  N  ( A  l 
CP=GOS(P.) 

SP*SIN(P> 

cf=rosio 

SfsSIM(C) 

C0*C0S I D) 

St)*SiN(D» 

CF*COS(F) 

SF«SIN(F) 

UC=SO“T(nSQfA*A )+A 
VO*DSO/UO 

ii]  =  sn»T(nso+C*C  Mf 
VI=DSQ/U1 

«J?=S0RT(nSQ4-R*P  )♦» 
V?sDSO/ll? 

U4=S0P  T( DCO+D»n)+D 
V4=D$Q/U4 

Uf*SQRT(DSQ«-F*F  J+F 
VODSQAJ6 

CALL  SICKST.CNUOI 

SU0*SI 

CU0=CI 

CALL  SIC  MS  I, Cl,  VO) 

SV0*SI 

CVO*CI 

CALL  SICICSI.CI.UU 

sm*si 

CU1*C! 

CALL  SICMSI,CI,V1I 

SVi’sS I 

CV1*CI 
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CALL  SICmi.CT.U2l 

$U2*SI 

ci»;*c  i 

CALL  SIC.ICSl.CI.V2) 

5V2*Sl 

CV2-CI 

CALL  SICKS  I.CI.U4) 

SH4*Sl 
CU4*G  I 

CALL  SICKSKCI.V4) 

$V4*S I 
CV4*C.I 

CALL  SIC  1 1S  T  ,CJ.,U6) 

su6*sf 

Cl!4=C  I 

CALL  SICnSI.CI.V4) 

SVfr*S! 

CV6*CI 

RL*DD4*CCA*(CUCtCV0n(S.U0-rSV0.L*SA).-D4*(CC*(CU!>CVl)>CR*au?'*-CV2K- 
TSC*I  SUl-$vi  nSP*C  SU2-SV2in>C6t<CU44GV4)>S0*ISU4-SVM>Ce*(:CUA+CV6) 
.♦♦SF*I'SUf>?$V6 ) 

AG=SE*(  CU4-CV6)— CF*(  Stl4*SV4)  +  S0*l  CU4— CV4I— Cn*rSU4+SV4)-P4*(  SC* 
*(CUl-CVl)>SB*ICU2-CV2)-CC*(SUi>SVn-C»»*fSiJ24SV2))-*-nn4*(SA*<CU0“ 
ICVO  j-CA*( SUO*SVO) I 
zvi-CMPtX(RL.Arr)*snK 
IFLN.60,.0).  RETURN 

A*v7N,tZ* 

R*A-rDK 

C*A+DK 

f\*P-r)K 

E=C +OK 
C  A*CO  S ( A  ) 

SA^SINt  A) 

CP*CO$ I B ) 

$P*$IN( P i 
CC=COS<C) 

SC*  SINK) 

cr=cos<ni 

SO*SINK) 

CE*COSTE) 

SE*S!N( F ) 

VO* SORT! OSO*P*D)*P 
UO*OSQ/VO 
Vl*SORTCOSO*B*BI+R 
Ul-OStWl 

V2*SQFT(DSO+A*A)+A 

U2*nSQ/V2 

V4*S0RTI0S0+C*C)*C 

U4*OSO/V4 

V6*S0RTCDSQ*E*F)<F 

U4*0S0/V6 

call  sicnsi.ci.uo) 

SUO*Sl 
CUO*C  I 


CALL  SICUST.C  I.VOI 

SVO»SI 

CVO*C! 

pin  ^ir.usi.ci*t»n 
$in*si 

Ct'l*CI 

CALL  SKMSUCI  .VI) 

SV 1 *S ! 

CV1«C  I 

(Ait  sicusi.r  i,u?i 
Sl'?=S  T 
eu2*c  i 

CALL  SICMSlt.CI.V?) 
sv?*s  l 
cv?*r.  j 

Call  5  IT  I  i  ST,  ci,  i  14) 

SU4=SI 

OI!4*C! 

CALI  SIC  ITS MCI ,V4l 
SV4=S  I 
cv4«r.  i 

CALL  S  IC  I‘(  S  T  ,C.I  *U6 ) 

SUf*St 
CU6=C  I 

CAIL  SIC I(SI »C! »V6 ) 

SV4=ST 
CVf.sC  I 

Ml*nn«*tCA*{.Ci|7«-CV?l-SA*!$U2-SV?)  )-D4*(CC*(  CH44-CV4)+C^*IC;1)I+CV1  )- 
«  ST*  (  SU4-SV4  )-SP*(  Sill-SVl ) )  +CO* ( CUO+CVO  J-SP*  I SUO-S VO )  +C E* ( CUA+C  V6 ) 
«-SF*LSHf~SVM 

AC*-SF*(  CIJA-C  VA)“CE*l  Sllfc+SVM-Sn*(TUO-rvO,)-Cn*I  SU0«-SV0K0.4*ISC*( 
.♦CU4-CV4  )  +  Srt*(  V.1M-CB*!  SllUSVl  U-CCX  SU44-SV4)  )>rin4*(-CA>M  Sl)?+ 

f SV? )-$A*( CH?-CV? )  ) 

7*r*cm  XHH,AG)*snK+7Ml 

R.f  TIJON 

ENP 


SHFCnilTIA'F  SICMSI,CI,X) 

C  »'"*  P I J  T  F  S  SINE  AND  Cf’STN'F  INTEGRALS 

7=ARS{ X) 
l  F  (  7—4,  )  1 , 1 , 4 
I  Y*(4.-7)*<4.«7) 

1  ST*X*(  (  Ml  1*7  S3141E“R*Y'*>  1,  *5f  H9R  RF-7  )  *Y+ 1,  T7  4i  AAE-S)  *v  +  6.<ne>RPQF-4) 
l*Y*l.RA4H«?F-M*Y+4,  JS,)SO<?F-l ) 

CI«((  S.77?15f.F-l+AlOG(Zn/7-7*(  MU  UTAAORSP-IO^YU,  5R49Qf>F-A)*Y 
1  +  1.  72*»7R?E-M*Y«T.  lASQQ9F-4)*Y>4,oqOl9?nF-A)  *Yfl.TlSRC«E-l )  )*7 
RETURN 
4  S!*SIN(?) 

Y  =  C.ns(7)  ' 

7.  *4  ,  /  7 

|)r  (  ( (  (  M  M4.04An49E-  4*7-?.?79l4TP-7)*7>S.5lS070E-?)*Z~7.?4U4?E-?) 


l*2>4.9877i6P-2»*Z-3.33?519E-3)*Z-2.3i46l7e-2>*2-i.i34<»5ftF-5)*2 

?♦>.  2*>001  IE-t2)*Z+7«  5R39FQF-10 

V#  C  <  C  ( I  (if  T>5.10R699K-3*?42.819i7<»F-2  l*7-6.53725XE->l  *7 
l>7.O02034F-2)*7-4.40r>4l6E-2)*Z-7.945556E-3>*7>?.60l293E-2>*7 
2-1 .  764000F-A )  *7-  3. i 2241 PF-2 1  *2-4. 644441  f- 7 ) *?.*? . 50000 JE-1 
Cl*?*  (SI*V-Y*IJ) 

Sr*-Z*CM*U4Y*V)  ♦1.570796 

«FT»PK' 

end 


SUPPPUTIMF  CSMIM  A,  N, DETERM) 

TNVFRTS  A  GIVEN  SQUARE  MATRIX 

CPMPt  EX  A ,  PI  VOT ,  AMA  X ,  T  ,  SWAP ,DF TERM, U, CMPl X , CONJ 

niMENSITN  I P! VOT ( 50) » INDEX! 50,  27 *  Al  50,501  ,PIVTT(50) 
OFT  ER  M*CMPL-Xf  1.0, 0.0) 
nr  ?o  j=i»n 

20  I P I  VO  T  (.J )  *0 

no  600  1*1,  N 
AMAX*CMPLX( 0.0,0.01 
no  105.  J*  1,  N 

IFMPIVnT(J)-l)  60,105,60 
60  no  100  K*i»N 

IM  IPIVnT(Kl-L)  80,100,740 
80  -TPMP*  AMA  X*OON  J  <  AMAX)-A<  J,K  )  *f.r)NJ  (AlJ,*)) 

1F(TEMP»85, 85,100 
85  TFOW=J 

IfriLIIMaK 
<\MAX*A(  J  ,K  1 
100  C  r  6*  T I  Nil  IF 
105  Cl  NT  INI  IF 

I  r  ?  vnT  (I.cm  UM » *  I  p  I  VOT  (I  COL  UM | ♦! 

IFUenu-lCOU/M)  140,260,140 
140  DETER  M*-HE TERM 
OP  700  1=1, M 
SWAP*AC IRHW.l  ) 

A(!ROW,L)=AI ICOLUM, L) 

200  AnC0l«»M,LI=SWAP 
260  IMngX< I,11=!R0W 

INDEX ( I , 2)*IC0LUM 
P IVPT  t I )  =  A( ICOLUM, ICOLUM) 

DETER  M*DFTEPM*PivnTU» 

TFMP«PIVnTm*CONJ  (PIVOTU)) 

IF(TEMP)  330,720,330 
’•<30  A(ICnu»M,  ICOLUM) *CMPLX<UO, 0.0) 
no  ^50  L*1 , N 

ll=P  IVOT  (  I  ) 

350  A< ICOLUM, l )*AlICnUJM,L)/U 
380  nr-  550  ti*l,N 

I F  < LI- ICOLUM 1 400, 550, 400 


51 


*C0  T*A(ll.!COLt!“l 

Atll,ICClUMI*CMPLXtO.C,0.0» 

00  450  t«l*N 

450  ACLi,l  )»A|L1.H-U*T 
55o  crMiwuF 
400  CCNTlNliF 

nc  710  I*  1^ 
t.=M  +  l-I 

i f ( iNofcxu. iHTNotxa.?) i  6*0,710. <*?o 

430  .iPr,wsf.wnFx(L,n 

jroLU^*lNn£X(L,?» 

DC  7C«  k«l,N 

SWAP=A(K,JRnW> 

A(k ,JHrw)=4(K* jCOtUMl 
A(Kf  JC.ni  I»W)*SWAP 
705  C  TNT  INI  10 
710  CONTINUE 
RFTU°M 

7?0  WP  I  T.F  (  3,730 1 

7*0  Fnf?MAT.(«C»,5X,»  MATR I X  IS  SINGULAR  •//) 

740  RETURN 
FNO 


FUNCTION  CPNJ(Z7) 

complex  r.‘'Pix,77,cnNj 

CPNJ«r.«PLXl  PEAl.(77),-AlMAGt77n 

RFTIIPN 

FNO 


Sl'PftnuTIMF  MIILTPY(L,M.N,A,R,C) 

C 

C  ‘MILT  |  PL  IF S  TWO  MATRfCFS 

C 

COMPLEX  A|50,50),R(50,l),CI50,ll 
nn  1  !  *  1 1 L 
no  IKs  1 »  N 
C ( I ,K)«C.C 
no  1  j=i,r 

1  CU,K|tf(  I,K  ) +A(  I  ,  J|*B{  J  ,KI 
RETURN 
END 
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NORMALIZED  F-  FIELD 

PATTERN 

0 

0.000001 

53.712891 

5 

C. 001 483 

66.059052 

6 

0.005793 

63.462265 

O 

0.012978 

59.139740 

1? 

0.022916 

53.071182 

15 

0.034086 

45.249466 

18 

6.048338 

35.656738 

21 

C. 061772 

24.258026 

24 

6.073137 

10.986127 

27 

0.079783 

-4.302690 

30 

0.078781 

-21.962402 

33 

0.067565 

“43.012985 

36 

0.045304 

“71.680725 

39 

0.02048? 

-141.776978 

42 

0.042617 

118.120651 

45 

0.083445 

80.080215 

48 

0.115672 

50.298401 

51 

0. 126?59 

21.C32181 

54 

0.105938 

-10.648551 

57 

0.053941 

-54.991577 

60 

0.044821 

161.263107 

6* 

0.129797 

106.062408 

66 

0.198734 

70.245300 

69 

0.21486? 

36.0554?0 

72 

0.151932 

-0.682838 

75 

0.030137 

-121.462845 

7* 

0.236566 

131.869217 

81 

0.503893 

95.018326 

84 

0.755990 

60.819748 

87 

0.935622 

77.297760 

90 

1.000000 

-5.863531 

93 

0.9334?! 

-38.657883 

96 

0.752121 

-70.907059 

99 

0.499261 

-102.105713 

102 

0.232193 

-130.074353 

105 

0.029792 

-81. 25«Oli 

108 

0.153752 

-29.985626 

111 

0.215196 

-55.546463 

114 

0.197901 

-82.309097 

117 

0.128347 

-105.952866 

120 

0.043571 

-107.648117 

123 

0.055012 

-20.748306 

126 

0.106574 

-31.308990 

129 

0.12636? 

-52.028732 

132 

0.115269 

-72.925568 

135 

0.082777 

-90.934418 

138 

0.041938 

-97.939575 

141 

0.020670 

-39.324295 

144 

0.045762 

-10.939811 

147 

0.067850 

-19.835052 

1*50 

0.078877 

-33.298767 

1*5  ^ 

0.079726 

-47. 126251 

156 

0.072073 

-60.234039 

159 

0.06159? 

-72.199876 

l**? 

0.048101 

-82.808686 

16K 

0.034663 

-91.933441 

168 

0.022626 

^99.493362 

171 

0.C12829 

—  j  06,  47-0761 

174 

0.005694 

“109. 682693 

177 

0.001 *00 

-112.246414 

190 

o.c jo:oi 

-96.120089 
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APPENDIX  C 

PROGRAM  FOR  WIRE  OVER  RADIAL  WIRE  SYSTEM 


This  program  is  suitable  £or  analyzing  radiation  by  wires  over  radial 
wire  systems  as  exemplified  by  Fig.  15.  llie  wires  are  assumed  lossless 
with  no  externally  applied  loading  and  excitations  are  restricted  to 
center-wire  positions  corresponding  to  peaks  of  triangle  current  expansion 
functions.  The  first  triangle  begins  at  the  junction  and  extends  towards 
the  end  of  the  center  wire.  The  junction  Itself  corresponds  to  the  peak 
of  the  first  triangle  of  the  branch  wire..  The  sample  output  is  for  a  X/4 
center  wire  with  two  branches,  each  of  length  X/3.  The  branches  are  normal 
to  the  center  wire  as  in  Fig.  15a.  The  center  wire  has  six  triangle  ex¬ 
pansion  functions  while  each  branch  supports  four.  The  wavelength  is  one 

-4 

meter  and  each  wire  is  of  radius  10  meter.  The  conflguraitcn  has  a 
single  unit  excitation  and  that  is  located  right  at  the  junction  Itself. 
Since  the  junction  marks  the  position  of  the  first  branch  triangle  which 
is  the  seventh  triangle  the  feed  point  is  assigned  the  number  seven.  Pat¬ 
terns  are  computed  for  the  $  *  0°,  90®  planes.  x,y,z-coordinates  of  axial 
points  defining  the  problem  geometry  are  labeled  PX,  PY,  PZ,  respectively 
in  the  sample  output. 

Subroutines  are  listed  first,  followed  by  the  main  program  and 
typical  output.  The  program  is  described  in  Chapter  4.  Further  informa¬ 
tion  on  computational  procedures  is  available  elsewhere  [4]. 


SUBROUTINE  UNEQ  <N,Z) 

COMPLEX  Z(  400 ) »  STOP ♦ STO  f  ST  t  S 
DIMENSION  LA ( 50 ) 

DO  20  I =1 »N 
L  A  (  I  )  =  I 
20  CONTINUE 
Ml  =  0 

DO  1R  M=l,N 
K=M 

DO  2  I =M,N 

K1=M1+I 

K2*M1+K 

IF  (CABSIZ(Kl)  )-CARS(Z(K2m  2»2tA 
A  K=  I 

2  CONTINUE 
LS=LA(M) 

LA(M)=LA(K) 

LA(K )*LS 
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K2=M1+K 
STflR=7  I K2 ) 

J1=0 

DO  7  J  =  1  ,N 
Kl= Jl+K 
K2=41*M 
ST0*ZIK1) 

Z(K1)=Z(K2) 

Z I K2 ) =STO/STOR 
Jl-Jl+N 

7  CONTINUE 
K1=M1+M 

Z ( K1 ) =1 «/STOR 
no  u  I* un. 

IFII-M)  12,  lit  12 
1?  K1*M1+I 
ST=Z(K1 ) 

Z(K1)=0. 

J1=0 

no  lo  4=i, n 
K2=J1+M 
Kl=4 1+1 

Z(K1)=Z(K1 ) — Z ( K  2 ) *ST 
J1=J1+N 

10  CONTINUE 

11  CONTINUE 
M1=M1+N 

1H  CONTINUE 
J1=0 

00  9  J  =  i,N 

IF  ( J-LA I 4  )  )  1 A , fi ,  14 
14  LAJ=l.AU) 

J2=(LAJ-1 )*N 
31  00  13  1  =  1, N 
K2= J2+I 
K1=J1+I 
S=Z(K2) 

Z(K2)*Z<K1) 

Z(K1)=S 
13  CONTINUE 

LAI J)*LA(LAJ) 

LA(LAJ)  =LAJ 

IF  I J-LAI 4 ) )  14,8,14 

8  JUJ1+N 

9  CONTINUE 
RETURN 
END 

SUBROUTINE  ROW! N, TH, PH, E ) 

COMPLEX  C 1 20 ) ,E I ? )  ,U ,U  1  »U2 ,U3 ,U4 ,U5 

COMMON  XX { 50 ) , X Y I  50 ) , XZ I  50 ) , TX I  50 ) , TY I  50 ) , TZ 1 50 ) , AL I  50 ) 
COMMON  T  ( 100) ,  TP  1 100)  »RAD2(  10),L(10),LL(10),LR(10),RK 
COMMON  /COA/  C 

OIMENSION  BKR I  50), DTI  50), DPI  50) 

U=(0.,i.) 

ETA*374,7307 
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CT=COS(TH) 

St=SIN( TH) 

CPsCOS(PH) 

SP=SIN(PH) 

S1=CT*CP 

S2=CT*SP 

BK1‘=BK*ST*CP 

BK2=BK*ST*SP 

RK3=BK*CT 

N2  =  l 

N3=-2 

DO  1  NS=l tN 
IF  <l(N2)-NS>  2  »3»  2 

3  KK=1 
N3=N3+2 
N2=N2+1 
GO  TO  4 

2  KK=3 

4  DO' 5  K=KK,4 
N7=N3+K 

BKR  ( N7 ) =XX ( N7 ) *RK 1+XY ( N7 ) *BK2+XZ ( N7 )  *BK3 
OT ( N7 ) =TX { N7 ) *S1 +TY ( N7 ) *S2-TZ ( N7 ) *ST 
DP ( N7 ) =-TX ( N7 ) *SP+TY ( N7 ) *CP 

5  CONTINUE 
N3=N3+2 

l  CONTINUE 
N2=  1 
N3=-2 
U3=0. 

04=0. 

DO  A  NS=ltN 
IF  (L(N2)-NS)  7» fl ♦ 7 
B  N3=N3+2 
N2=N2+l 
7  JI=(NS-1I*4 
01=0. 

02  =  0. 

DO  9  JS*lt4 
J2= J1 +JS 
J3=N3+.IS 
S 1  =  BKR  (  J  3 ) 

05  =  T  ( .12 )  *  (COS  (  SI  )+0#SIN(  SI ) ) 

Ul  =  Ul+Ui>*DT(J3) 

U2=U2+U5*DP( J3) 

9  CONTINUE 

03=U3+01*C ( NS ) 

U4=U4+U2*C(NS) 

N3=N3+2 
A  CONTINUE 

S1=.0795774*ETA*BK 

E( l )=-U#Sl*U3 

E(2)=-U*S1*U4 

RETURN 

END 
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SUBROUTINE  GOMTR YINBRCH, NEC ,NEB»N1»N»CLENTH»BLENTH» THETA) 

THIS  SUBROUTINE  GENERATING  PARAMETERS  OF  A  UMBRELLA  SHAPE  WIRE 

STRUCTURE  (ANTENNA  OR  SCATTEREE) 

I/P  NBRCH=  NO.  OF  BRANCHES 

NlsTOTAL  NO.  OF  SEGMENTS 
N= TOTAL  NUMBER  OF  EXPANSION  FUNCTIONS 
NEC=NO.  OF  EXPANSION  FUNCTION  OVER  CENTRAL  WIRE 
NEB*NO.  OF  EXPANSION  FUNCTION  OVER  EACH  BRANCH 
CLENTH=LENGTH  OF  CENTRAL  WIRE  MAY  BE  NEGATIVE  TO  REPRE¬ 
SENT  THE  CASE  WITH  NEGATIVE-DIRECTED  CENTRAL  WIRE 
BL6NGTH=LENGTH  OF  EACH  BRANCH  UN  METER) 

THETA=ANGLE  BETWEEN  EACH  RR  ANCH  AND  THE  POSITIVE  Z^AXIS 

fl/P  XX(I)  =  PTH  COMPONENTS  OF  THE  COORDINATES  OF  THE  CENTER 
OF  THE  I'TH  SEGMENT 
XY (  I )  I  S  THE  Y  COMPONENTS 
XZ  IS  THE  Z  COMPONTS 
AL(I)  IS  THE  I'TH  SEGMENT  LENGTH 
TX  IS  THE  PROJECTION  OF  AL  TO  X-AXIS 
TY  IS  THE  Y  COMPONTS 
TZ  IS  THE  Z  CONPONENT 

COMMON  XX ( 50) »  XY ( 50 ) , XZ ( 50) ,TX ( 50 ) ♦ TY ( 50 ) , T7 ( 50 ) , AL ( 50 ) 

COMMON  T ( 100) *TP ( 100) ,RAD2< 10) ,L ( 10 ) ,LL ( 10 ) ,LR ( 10) 

DIMENSION  PX{ 50 ) , PY ( 50 ) »PZ( 50) 

DOUBLE  PRECISION  DP1,0P2,DP3 

NPR=2*NEB+3 

NPC=2*NEC+3 

DP1=NPC-1 

DP2=CLENTH 

PP2=DP2/DP1 

L  ( 1 )  =1 

LL ( 1 ) =1 

1(2) -1+NEC 

LL ( 2 ) =1+NPC 

DO  10  1*1, NPC 

Si=I~l 

PX(I ) =0 . 

PY(I )=0. 

PZ (I )=Sl*DP2 
10  CONTINUE 
DP1=THETA 
DP2-.G174533 
TH=DP1*DP2 
DP1=BLENTH 
DP3*2*NEB 

DP3=DP1/DP3 

D3=DP3*C0S(TH) 

D1=DP3*SIN(TH) 

Jl=NPC+l 
DO  20  I =1»NBRCH 
RRCH=NBRCH 
X=  I -1 
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0P1=3A0.*X/BRCH 

PH=0P1*DP2 

02=01*SIN(PH) 

h4=ni*cns(PH) 

LL  C  l +2 ) =LL ( I +1 ) +NPH 
L  (  I+2I=LU+1)+NEB 
no  40  .1  =  1*3 
PX(  Jl)=0,. 

P V ( J 1 ) =0 • 

PZ(  J1  )  =PZ(4-.I) 

Jl=Jl+l 
40  CONTINUE 
NP1=2*NER 
00  30  J=l«NPl 
S1  =  J 

PX( J1 )=04*S1 
PY(Jl)=n2*Sl 
P Z  <  J 1  )  =  03*51 
J1=J1+1 
30  CONTINUE 
20  CONTINUE 

L ( MBRCH+2 ) =200 

LL ( NBRCH+2 ) =200 

NP=NPC+NBRCH*NPB 

WRITE  (  3*301 )  (PXU  )*I=1»NP) 

WRITE  (3r302)  < PY ( I ) * I = 1 t NP ) 

WRITE  (3*303 ) ( PZ ( I) t I =1 tNP ) 

301  FORMAT (  ' OPX  •  / ( IX*  7EU.4H 

302  FORMAT ( • OPY • / ( IX  *  7E11.4)) 

303  FORMAT! 'OPZ '/(IX*  7EU.4)) 

.11=1 

J2  =  2 
N1=0 

DO  2  J  =  1 » NP 
IF  (LL(Jl)-J)  3,4*3 
4  J?= J2-1 
L( J1)*J2 
J1=J1+1 
GO  TO  2 
3  N1=N1+1 
J3=J-1 

IF((Nl/2*2-Nl).EO.O)  J2  =  J2-t  ’ 
XX(N1)=*5*(PX(JI+PX(J3)> 

XY(N1 )  =  . 5^ ( PY  < J)+PY( J3 ) ) 
X7.(N1)  =  *5*(PZ  ( J)  +PZ  (  J3 ) ) 
S1=PX(J)-PX(J3) 

S2=PY(J)-PY( J3) 

S3=PZ ( J  )-PZ ( J3 ) 

S4=S0RT(S1*S1+S2*S2+S3*S3) 

TX(N1)=SI/S4 

TY(N1)=S2/S4 

TZ ( Nl )=S3/S4 

AL ( N 1 )=S4 


?  CONTINUE 
L ( J1 ) *J2 
N»J2-2 
Jl*l 
J2*-2 

on  11  J=1  ,N 

IFM(Jl)-J)  13 » 14  ♦IB 
14  J2*J2+2 
J1*J1+1 
13  J3=(J-1)*4 
J4=J3+j 
J5-J4+1 
Jft* J5+1 
J7=J6+1 
K4sJ2+l 
K5=K4+1 
K(S=K5+1 
K7=KA+1 

S1=AL (K4)+AL(K5 ) 

S2sAL ( Kft ) +AL ( K7 ) 

T  <  J4V1*  AL  ( K  4 )  *  .  5*AL  (K4)  /SI 
T( J5)=AL(K5)*{AL(K4)+.5*AL<K5) )/Sl 
T(  J6)=AL(KA)*l  AUK7)  +  .5*AUKM  )/S2 
T( J7)=ALiK7)*.5*AL(K7)/S2 
TP( J4)=AL(K4)/S1 
TPIJ5)=AL(K5)/S1 
TP< JA)=-AL(KA)/S2 
TP(J7)=-AL(K7)/S2 
J2= J2+2 
11  CONTINUE 
RETURN 
END 


SUBROUTINE  Z RRCH (NBRCH »N1,N,NEC»NEB»Y) 

THIS  ROUTINE  COMPUTE  THE  GENERALIZED  IMPEDANCE  MATRIX  OF  UMRRELLA 

SHAPE  OF  WIRE  STRUCTURE 
• 

I/P  NRRCH=  NO.  OF  BRANCHES 

Nl=TOTAL  NO.  OF  SEGMENTS 

N=TOTAL  Nl MRER  OF  EXPANSION  FUNCTIONS 

NEC=NO.  OF  EXPANSION  FUNCTION  OVER  CENTRAL  WIRE 

MFR=MO.  OF  EXPANSION  FUNCTION  OVER  EACH  BRANCH 

XX,XY,XZ  ,AL»TX,TY,TZ  ARE  GEOMETRY  PARAMETERS 

RAO=RADIUS  OF  EACH  WIRE 

O/P  Y  IS  THE  GENERALIZED  IMPEDANCE  MATRIX  COMPUTED 

COMPLEX  Z(  400), Y<  400),  PSI  (  200  )  ,IJ,U1  ,U2,U3,U4,IJ5,UA 
COMMON  XX ( 50 ) »XY ( 50) , XZ < 50 ) ,TX I  50 ) ,TY ( 50 ) , TZ ( 50) , AL ( 50) 

COMMON  TUOO),TP<  100 )  ,R  AD2(  10)  ,L  ( 10) ,  LL  ( 10)  ,  LR  ( 10 ) ,  RK 
DIMENSION  DC  I  200 ) 

U=(0.,1.) 

P I =3.141593 
FTA=37A.7307 


62 


Cl* • 125/PI 

C2*.25/?I 

U3«U*RK*ETA 

U4*-U/BK*ETA 

BK2»RK*RK/2. 

RK3=BK2*BK/3. 

NP*0 

M2*l 

NO=l 

N3=-2 

mUnec+nfb 

M2*2*Ml*4 
M4*NEC 
M3*2*NE02 
no  10  NS*  1  *H,1 

I F(NS.GT.NEC)  M3=Nl-( (NBRCH-1 ) /2)*< 2*NEB+2 1 
IF(NS.GT.NEC)  M4*NEC+< (NBRCH+2 ) /2 ) *NEB 
IF ( L ( N? )-NS )  12 1 1 1 *12 
11  KK*1 
N3=N3+2 
N2=N2+1 
GO  TO  13 
1?  KK*3 

00  14  NF=1,M3 

N4=NF+M3 

N5*N4+M3 

NA=N5+M3 

DC!  NF  )  =  f)C(N5) 

DC.(N4)=DC(NA) 

PSl (NF )=PSI (N5) 

PSI (N4)=PSI (NA) 

14  CONTINUE 

13  IF  ( N3+1-LR ( NO) )  5»At5 
A  AA=RAD2 ( MO ) 

N0=N0+1 
5  CONTINUE 
DO  15  K=KK*4 
N7=N3+K 
Kl* C K-l  )*.“? 

DO  1A  NF=1 »M3 
Nft=NF+Kl 

S1=XX ( N7 )-XX ( NF ) 

S?*XY(N7 )-XY (NF) 

S3=XZ(N7)-XZ(NF> 

R2=S1*S1+S2*S2+S3*S3+AA 

R*S0RT(R2) 

RT=ABS(S1*TX(N7)+S2*TY(N7)+S3*IZ(N7) ) 

RT?*RT*RT 

RH= ( R2-RT2 ) 

ALP*. 5*AL ( N7 ) 

AR=ALP/R 

S1=BK*R 

II2=C0S(S1  )-*U#SIN  (SI) 

IF(AR-.l)  22*  22*21 
21  U2=U2*C1/ALP 
S1=RT-ALP 


- ■  iiiiiiMiiiMi . Mi  — .1.1. <■<.  *  .  l 
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S2=RT+ALP 

S3*S0R T  ( S I  *S,l  +RH ) 

S4*SQRT(S?*S2+RH) 

IF (SI)  18,18,19 

1R  A'l  l*ALDfi(  (S2+S4)*(-S1>S3)/RH) 
GO  TO  20 

14  AI  l*AL0GnS2+S4)  /( S1+S3) ) 

20  Af2*AL(N7) 

A I  3*  ( S2.4S4-S  l*S3+RH#AI  1 )  /  2  . 


AU»AI2*(RH+ALP*ALP/3.  +  RT2) 

S3=AI1*R 

S1*AI1-BK2*(AI3-R*(2.*AI2-S3)) 

S?=-BK*<  AI?-S3)+BK3'MAI4-3.*AI3#R+R2iM3.*AI2-*S3) ) 
GO  TO  28 
22  U?=U2*C?/R 
RA*BK#ALP 
BA2=RA*RA 
AR2=AR*AR 
AR3=AR2*AR 
7.R-RT /R 
ZR2=7R«7.R 
ZR3=ZR2*ZR 
ZR4=ZR3*7.R 

Hl=(3.-30.*ZR2+35.*ZR4)#AR3/40. 

Al=AR*(-l.+3.*ZR2)/A.  +  (3.-3n.’S‘ZR2+35.*ZR4)*AR3/4n. 
A0= 1 • +AR *A 1 

A2=-ZR2/A.-AR2*( l.-l?.*ZR2+15.*ZR4)/40. 
A3=AR*(3.*ZR?-5.*ZR4)M0. 

A4=ZR4/120. 

S1=A0+BA2*( A2+BA2*A4) 

S2=RA*(A1+RA2*A3 ) 

28  PSI  (NR)*U2*(Sl+ll*S2) 

DC (N8 ) =TX ( NF ) #TX (N7 ) +TY ( NF ) ^TY ( N7 ) +TZ (NF ) fcTZ ( N7 ) 
l  A  CONTINUE 
15  CONTINUE 
N3=N3+2 
J3= ( NS-1 ) *4 
J7  =  -2 
J9  =  l 

no  25  NF= 1 ,  M4 
Jl=(NF~l »*4 
IF ( L ( J9 )-NF )  2A,27,2A 
27  J9=J9+1 
J7=J7+2 
2 A  N9=N9+1 
U5  =  0. 

UA=0  • 

J5  =  0 

DO  23  JS= 1 ,4 
J4= J3+JS 
J  8  J  5  +  J  7 


DO  24  JF  =  l ,  4 
JA=JR+JF 
J2= Jl+JF 

U5  =  T(  J2)*T(J4)*DCUA)*PSI  ( JA)+IJ5 
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I  IA= TP  1 .1? )  *TP  ( J4  )  *PS  I  UA ) +UA 

24  -CONTINUE 
J5=J5+M3 

23  CONTINUE 

7 IN9)=U5*03+UA*U4 
,17  =  J7+2 

25  CONTINUE 
10  CONTINUE 

J1  =  0 
J2=0 

M3=NEC*NFC 
no  30  1=1 »NEC 
,13*M3+I 
00  31  J=1 »NEC 

.12  =  J2+1 
.ll=Jl  +  l 

yi  ji)=zur.) 

31  CONTINUE 
S1=NR»CH 

no  3?  .1  =  1 1 N  E  B 
Jl=Jl+l 

Y{ J1 )=Z ( J3)*Sl 
J3= J3+M4 

32  CONTINUE 
30  CONTINUE 

J2=M3 

no  35  I =1 » NEB 
no  3A  J=1,NEC 
.12  =  J2  +  1 
J 1 = J 1+1 
Y ( J 1 ) =7  C  J?  ) 

3A  CONTINUE 

00  37  ,1  =  1  »NFR 
Jl=Jl+l 
J2=J2+1 
Y  C  J 1 )  =  7  (  J  2 ) 

I F ( NBRCH , EO • 1 )  GO  TO  37 
KK=(NRRCM+?) /? 

DO  41  K  =  2 » KK 
KKK= ( NRRCH/2 ) *2-NBRCH 
I F ( (K.FQ.KK) .AND. (O.EO.KKK) )  GO  TO  42 
J3=J2-MK-1)*NFR 
Y( J 1 ) =Y ( J1 )+7 ( J3>*2 
00  TO  41 

4?  J3=J2+(K-1 )#NER 
Y(J1)=Y( J1)+Z(J3) 

41  CONTINUE 
37  CONTINUE 

J?=J2+(NBRCH/?)*NEB 
35  CONTINUE 
RETURN 
ENO 
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MAIN  PROGRAM  MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM 

COMPLEX  7  (.400)  »U(  20)  ,C( 20) «E(3)  »EI ( 2) ,UV (4) »U1,V 
COMMON  X  X ( 50  > ♦ X Y ( 50 ) , X  Z ( 50 ) , T  X ( 50  > ,T Y  < 50  j  ,  T7  ( 50 )  ,  AL  (  50 ) 

COMMON  T(100),TP( 100) ,RAD2( 10) ,L< 10) ,LL( 10) ,LR< 10) ,BK 

COMMON  /COA/  C 

0!  MENS!  ON  I FP ( 20) ,RAD( 10) 

1  READ( i » 101 , END* 5 00)  NB,NFB,NEC ,BK,BLENTH,CLENTH, ANGLE 
101  F0RMAT(3I3,3E14.7,F6.2) 

WRITE  (  3»301:)NB,NEB»NEC,RK 

301  FORMAT ( '0  NB  NEB  NEC  BK* /IX , 13 » 2 14, E14.7 ) 

WRITE  ( 3 f 302 )  BLENTH»CLENTH, ANGLE 

302  FORMAT!  'OLENGTH  OF  EACH  BRANCH  »/lX,E14.7//  •  LENGTH  OF  CENTER  WI 
1RE'/1X,F14.7//'  ANGLE* /IX  ,F7»2) 

NR  =  1 

LR ( 1 ) =1 

LR (NR+1 )=200 

RE AD (1,109)  RAO! 1 ) 

109  FORMAT (E14.7) 

WRITE  (3,305)  R AD ( 1 ) 

305  FORMAT! *0RAD'/1X,E14. 7) 
no  46  1=1, NR 

RAD2!  I  )=RA()(I  )*#2 
46  CONTINUE 

CALL  GOMTRY ( NR » NEC, NER, N1 ,N , CLENTH, BLENTH,  ANGLE ) 

CALL  ZBRCH(NR,N1,N,NEC,NEB,Z) 

NE=NFC+NFR 

NE2=NE*NE 

CALL  LINFQ  (NE,Z> 

DO  45  1=1, N 
U( I )=0. 

45  COMTINUF 

READ  (1,107)  NF 
107  FORMAT (2013) 

WRITE  ( 3,303)NF 

303  FORMAT  ( • ONF * / 1 3 / • OFP  VOLTAGE') 

DO  44  I = 1 , NF 

READ  (1,119)  Ji,V 

119  FORMAT! I3,2E14.7) 

WRITE  (3,120)  J1,V 

120  F0PMAT(IX,I3,2E14.7) 
l)(  J1  )=V 

I FP ( I  )=J1 
44  CONTINUE 

WRITE  (3,111) 

111  FORMAT  ( 'OCURRENT'/'  1  REAL  IMAG  MAGNITUDE 

1  PHASE') 

J3=0 

J1=0 

5  J1=J1+1 
1)1*0, 

J2=0 

6  J2= J2+1 
J3=J3+) 

Ul=Ul+Z( J3)*U( J2> 

IF  (J2-NE)6,7,7 
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nn  <»  i*i*ne 
m«c<n 

CM*CARS(U1) 
tF  (CM)  11  * i?« 1 1 
n  CP*A7AN2CA1MAG(U1),RFAL(UI))*57,2«5R 
GO  TO  10 
l?  CP=0. 

10  CONTINUE 

WRITF  (  3 1-1 1 2 )  1  tC(I )  *CMtCP 

112  FORMAT  (I5»3E14.6*F10.3) 

Q  CONTINUE 

I F (NR • FO  « 1  )  GO  TO  A? 

Jl=NER+NFC 
OH  41  I =2* NR 
J2=NER 

00  41  J=ltN6R 
J1=J1+1 
J2= J2+1 
C(J1)=C(J2) 

41  CONTINUE 

42  CONTINUE 

WR I TF  ( 3 « 1 13 ) 

113  FORMAT  ( *  OF 1 ELO  PATTERN*/'  K  =  1  FOR  THETA  COMPONENT  K  =  2  FOR  PHI 


1  COMPONENT' /'O  THETA  PHI  K  REAL  IMAG  MAGNITUD 

2E  PHASE') 

no  17  IPH=1,2 
PHIM  IPH-DwPO 
PH=PHI*. 0174533 
00  17  I TH  =  1  *  18 1 »  5 
THF=ITH-1 
TH=THE*. 0174533 
CALL  ROW  (N»TH,PH,E) 

00  16  K=1 »? 

Ul=E(K) 

GM=CABS(U1) 

IF  (GM) 13» 13» 14 

14  EP=ATAN2  (AIM  AG  (111)  *  RE  AL ( H 1 )  ) 

GO  TO  15 

13  EP=0 • 

15  GP=EP*57.2A5R 

WRITE  (3*114)  THEtPHl tK,Ui,GMtGP 
114  FORMAT (F6.0,F5. Ot 1 3 ♦ 3fi 14. 7 , F 10 . 3 ) 

16  CONTINUE 

17  CONTINUE 
GO  TO  1 

500  STOP 
END 

/£ 

//GD.SYSIN  nn  v 

2  4  6+0.6283185E+01+0. 1250000E+00+0. 2500000E+00+90. 00 

+0. 1000000 E-03 
1 

7+0.1000000F+01+n.0000000E+00 
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NB  NFB  NFC  BK 
2  4  6  0.6283184E  01 

LFNGTH  OF  EACH  BRANCH 
0.1250000E  00 

LENGTH  OF  CENTER  WIRE 
0.2500000F  00 


ANGLE 

40.00 


RAD 

0.9999999E-04 

PX 


o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 1563E-01 

0.3125E-01 

0.46R8E 

0.6250E-01  0.7813E-01  0.9375E-01  0.1094E  00  0.1250E  00  0.0  0.0 

0.0  -0.l563E-01-0.312f>F-01-0.46RRE-01^0.6250E-01-0.7813E-01-0.°375F-01 

—0.-104 4E  00-0.1250E  00 


PY 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.50R4E-08 

-0.1018E 

-07-0. 1527E-07 

-0.2036E 

-07-0. 2544E-07-0. 3053E 

-0.3542E-07-0.4071E-07 


PZ 

0.0  0.17B6F-01  0.357IE-01  0.5357E-01  0.7I43E-01  0.B929E-01  0.1071E  00 

0.1250E  00  0.1424E  00  0.1607E  00  0.17R6F  00  0.1964E  00  0.2143E  00  0.2321E  00 

0.2500E  00  0.3571E-01  0.1786E-01  0.0  O.A906E-0*  0.9812F-0R  0.1472E-07 

0.196PF-07  0.2453E-07  0.2944E-07  0.3434F-07  0.3925E-07  0.3571E-01  0.17R6E-01 

0.0  0.490AF.-0R  0.4812E-08  0.1472E-07  0.1962E-07  0.2453E-07  0.2944E-07 

0.3434E-C7  0.3925E-Q7 

NE 

1 

FP  VOLTAGE 

7  0 *  1 OOOOOOF  01  0.0 

CURRENT 

i  REAL  imag  magnitude  phase 

1  -0.51 1069F-03  -0.476440E-02  0.479223E-02  -94.105 

2  -0.50R337E-03  -0.437440E-02  0.4403R4E-02  -94.612 

3  -0.471473F-03  -0 .380«23E-0?  0.3R3R29E-02  -97.039 

4  -0.400856E-03  -0. 307407E-02  0.310009E-0?  -97.412 

5  -0.29R497E-03  -0.21RR66E-02  0.220R92E-02  -97.749 

6  -0. 167006F-03  -0. 117641E-02  0.11RR20E-02  -9R.063 

7  0.241R28E-03  0.253185E-02  0.254337E-02  84.529 

R  0.195378E-03  0.19414RE-02  0.200105E-02  84.382 


9  0.138629F-03  0.139732E-0?  0.14041RE-02 

1  A  A  *9  1  r  1  A  y  W  A  J  A  i  J  A  /  A  r  A  ft  A  ■»  i  A  f*  A  «  f  A  A 


84.319 


FIELD  PATTERN 


K= 1  FOR  THFTA  COMPONFNT  K=2  FOR  PHC  COMPONENT 

THFTA  PH!  K  REAL  I  MAC,  MAGNITlJE  PHASE 

0.  0.  1  0n852448F-01-0.2404486E-02  0. 18679R7F-01  -7=394 

0.  0.  2-0.43 1 84 37 E -OR  0 .23 1 87.54R-09  0.432465RE-0R  176. R96 

5.  0.  1  0*?fl 1 1 516E-01  0.4627321F-02  0.2R49341E-0]  9.345 

5.  0.  2-0.4283006E-08  0.541097RF-09  0.431 T048E-0R  172.770 

10.  0.  1  0.3760987E-01  0. 1 145996E-01  0. 3931 709F-01  16.943 

10.  0.  2 -0 .423 279 2E -OR  0.8457961E-09  0.4316465E-0R  168.671 

15.  0.  1  0.4702598E-01  0. 17R809RE-01  0. 5031075E-01  20.R15 

15.  0.  2 -0 • 4 1688 04 E -OR  0.1143O16E-0R  0.4322661E-0R  164.439 

20.  0.  1  0.5636786E-01  0. 236R479E-01  0.6114167F-01  22.787 

20.  0.  2-0.4092520E-08  0. 1429995E-08  0.4335156E-0R  160.712 

24.  0.  1  0.65621 2 6E -01  0. 2867664F-01  0.7161349t-01  23.601 

25.  0.  2-0 .4005798 E -OR  0.1704222F-08  0.4353250E-0R  1.56.926 

*0.  0.  1  0.7474923E-01  0. 3267739E-01  0.8157974E-01  23.609 

30.  0.  2-0.39 l 0849 E-08  0.1963477E-0R  0.4376066F-0R  153.314 

35.  0.  I  0.8368844E-01  0. 3552905E-01  0.90917R9F-01  22.999 

35.  0.  2-0.3R  101 15 E-08  0.2205868E-0R  0.4402590E?-OR  149.905 

40.  0.  1  0.9234822F-01  0.3710096E-01  0.9952223E-01  21.8R4 

40.  0.  2-0.37062 18E-08  0.2429835F-0R  0.4431719E-0R  146.725 

45.  0.  1  0.1006104E  00  0. 3729643E-01  0.i07300RE  00  20.336 

45.  0.  2-0.36018 69 E-08  0.2634155E-0R  0.446231 1E-0R  143.796 

50.  0.  1  0. 1083333E  00  0. 3605935E-01  0.1141769E  00  18.407 

50.  0.  2-0.349978 IE-08  0.2817917E-0R  0.4493231E-0R  141.135 

55.  0.  1  0. 1153556E  00  0. 3338032E-01  0.1200R80E  00  16.136 

55.  0.  2-0.3402594E-08  0.29R0492E-08  0.4523379E-08  13R.759 

60.  0.  1  0. 1215046E  00  0. 2930132E-01  0.1249R77E  00  13.556 

60.  0.  2 -0 .331279 0E-08  0.3121491E-0B  0.4551733E-08  136.679 

65.  0.  1  0. 1266061E  00  0. 2391 834E-01  0.128R456E  00  10.696 

65.  0.  2-0.3232655E-08  0.3240725E-0R  0.4577373E-08  134.905 

70.  0.  1  0.1304952E  00  0. 17 381R0E-01  0.1316477E  00  7.586 

70.  0.  2-0.3 1642 OOE -08  0.3338 133E-08  0.45994R5E-08  133.444 

75.  0.  1  0. 1330277E  00  0.9892R62E-02  0.1333951E  00  4.252 

75.  0.  2-0.3 1091 31 E-08  0.3413756E-0R  0.4617405E-08  132.303 

HO.  0.  1  0. 1340917E  00  0. 16975RRE-02  0.1341023E  00  0.725 

80.  0.  2-0.306R799E-08  0.3467675E-08  0.4630579E-0R  131. 4R5 

85.  0.  1  0.1336166E  00-0. 6923538E-02  0.1337957E  00  -2.966 

85.  0.  2 -0.3 044 19 8 E-08  0.34999R5E-08  0.463R643E-0R  130.993 

90.  0.  1  0.1315810E  00-0. 1566721E-01  0.1325104E  00  -6.789 

90.  0.  2-0.30359 32E -OR  0.3510747E-08  0.4641358E-0R  130. R29 

95.  0.  1  0.1280149E  00-0. 2422226E-01  0.1302B62E  00  -10.713 

95.  0.  2-0.3044 198 E-08  0.3499985E-0R  0.4638643E-0R  130.993 

100.  0.  1  0. 1230005E  00-0. 3228397E-01  0.1271666E  00  -14.704 

100.  0.  2 -0.3 068 7 99 E-08  0.3467677E-0R  0.4630582E-0R  131. 4R5 

105.  0.  1  0 . 1 166678E  00-0. 3956900E-01  0.1231953E  00  -18.732 

105.  0.  2-0.3 109 131E-0R  0.3413756E-0R  0.4617405E-0R  132.303 

ilO.  0.  1  0.1091857E  00-0.4582763E-01  0.1184132E  00  -22.765 

110.  0.  2 -0.3 164 2 OOE -08  0.3338 133E-0R  0.45994R5E-0R  133.444 

115,  0.  1  0.1007528E  00-0. 5085400E-01  0.112B594E  00  -26.777 

115.  0.  2-0 .323 26 57 E-08  0.3240725E-08  0.4577373E-08  134.905 

120.  0.  1  0.915R385E-01-0.5449336E-01  0.106569RE  00  -30.748 

120.  0.  2-0.33 1279 OF -08  0.3121491F-0R  0.4551733E-0R  136.679 

125.  0.  1  0.8189815E-01-0.566451 2F-01  0.9957R98E-01  -34.664 

125.  0.  2-0.3402594E-08  0.2980492E-0R  0.4523379E-08  13R.759 
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170.  0.  1-0. 273745RE-02^0. 15483686-01  0.15723806-01  -100.009 

170.  0.  2-0,42327926-08  0.84579706-09  0.4316465E-08  168.671 

178.  0.  1-0. 1072789E-01— 0. 66317366-02  C. 12612206-01  -148.251 

175.  0.  2-0.42R3006E-08  0.5M0976E-09  0.4317048E-08  172.770 

180.  0.  1-0, 18524506-01  0.24045146-02  0. 1R679R9F-01  172.574 

180.  0.  2-0. 43 184 376-08  0.23187626-09  0.43246586-08  176.896 

0.  90.  1  0.1497R44E-08-0. 52308316-09  0. 15R6554E-0R  -19.247 

0.  90.  2-0. 18524486-01  0.2*044866^0?  0..1R679R7E-01  172.574 

5.  90.  1  0.94915816-02  0.53235146-02  0. 10RR255F-01  29.282 

5.  90.  2-0.18524486-01  0.24044866-02  0. 1867987E-01  172.574 

10.  90.  1  0.19013026-01  0.10486356-01  0.21713096-01  28.873 

10.  90.  2-0.1852448E-01  0.24044R56-0?  0. 18679R7E-01  172.574 

15.  90.  1  0.285R743F-01  0.15328946-01  0.32437916-01  28.196 

15.  90.  2-0.18524486-01  0. 24044826-02  0.18679876-01  172.574 

20.  90.  1  0.38224806-01  0.19694566-01  0.4300012F-01  27.254 

20.  90.  2-0.1852448E-01  0.2404482E-02  0. 18679876-01  172.574 

25.  90.  1  0.47916016-01  0.23431566-01  0.53338376-01  26.055 

25.  90.  2-0. 18 52448 E-01  0.240448  IE-02  0. 18679R7E-01  172.574 

30.  90.  1  0.57628536-01  0. 2639677E-01  0. 63386386-01  24.606 

30.  90.  2-0. 18 52448 E-01  0.24044R IE-02  0.18679876-01  172.574 

35.  90.  1  0.6730169E-01  0. 2845993F-01  0. 73071 7RE-01  22.918 

35.  90.  2-0. 18 52448 E-01  0.2*04480E-02  0. 18679R7E-01  172.574 

40.  90.  1  0.7684582E-01  0. 2950929F-01  0.8231 694F-01  21.003 

40.  90.  2-0, 1852448 E-01  0.24044R0E-02  0.1867987E-01  172.574 

45.  90.  1  0.86141476-01  0.29457926-01  0.91039066-01  18.876 

45.  90.  2-0.18524486-01  0.2404479E-0?  0.1R679R7E-01  172.574 

50.  90.  1  0.9504282E-01  0.28250446-01  0.99152516-01  16.551 

50.  90.  2-0.185 2448E-01  0.2404479E-02  0. 1R679B7E-01  172.574 

55.  90.  1  0.1033820E  00  0. 2586960F-01  0.1065695E  00  14.046 

55.  90.  2-0.18524486-01  0.24044766-02  0. 1R679RRE-01  172.574 

60.  90.  1  0.1109769E  00  0. 2234155F-01  0.11320346  00  11.380 

60.  90.  2-0. 1852448 E-01  0.2404476E-02  0. 186798RE-01  172.574 


150.  90.  I  0.493831 7E -01-0. 388 44 R2E-01  0. 6283003F-01  -38.182 
150.  90.  2-0. 185 24 ARE -01  0.24044806-02  0. 18679R8E-01  172.574 
155.  90.  1  0. 40698226-01-0. 3371320F-01  0.52848126-01  -39.630 
155.  90.  2-0.1R5244BE-01  0.2404480E-02  0. IR679R8E-01  172.574 
160.  90.  1  0. 32222386-01-0. 2784979F-01  0.42589816-01  -40.R30 
160.  90.  2-0.18524486-01  0 .24044806-02  0. 1R67988E-01  172.574 
165.  90.  1  0.23952306-01-0. 2139971F-01  0.32119476-01  -41.771 
165.  90.  2-0.1P52448E-01  0,24044816-02  0. 18679886-01  172.574 
170.  90,  l  0.1585945E-01-0. 14510086-01  0.21495696-01  -42.449 
170.  90.  2-0.18524486-01  0.240448 16-02  0. 1R67988E-01  172.574 
175.  90.  T 0.7B95697E-02-0. 73279486-02  0. 1077222E-01  -42.857 
175.  90.  2-0.1852448E-01  0 ,240448 1 E-02  0. 1R&79RBE-01  172.574 
180,  90.  1-0.3092213E-07  0. 279622RE-07  0.416901-+6-07  137.854 
180.  90.  2-0.1852448E-01  0 .24044846-02  0. 1R679RR6-01  172.574 


